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Abstract 

Gershtein S.S., Kiselev V.V., Likhoded A.K., Tkabladze A.V. Physics of Be mesons: IHEP 
Preprint 94-81. - Protvino, 1994. - p. 82, figs. 23, tables 28. 

In the framework of potential models for heavy quarkonium the mass spectrum for the system 
(be) is considered. Spin-dependent splittings, taking into account a change of a constant for 
effective coulomb interaction between the quarks, and widths of radiative transitions between 
the (be) levels are calculated. In the framework of QCD sum rules, masses of the lightest 
vector B* and pseudoscalar Be states are estimated, scaling relation for leptonic constants of 
heavy quarkonia is derived, and the leptonic constant /s^ is evaluated. The Be decays are 
considered in the framework of both the potential models and the QCD sum rules, where the 
significance of Coulomb-like corrections is shown. The relations, following from the approximate 
spin symmetry for the heavy quarks in the heavy quarkonium, are analysed for the form factors 
of the semileptonic weak exclusive decays of Be- The Be lifetime is evaluated with the account 
of the corrections to the spectator mechanism of the decay, because of the quark binding into 
the meson. The total and differential cross sections of the Be production in different interactions 
are calculated. The analytic expressions for the fragmentational production cross sections of Be 
are derived. The possibility of the practical Be search in the current and planning experiments 
at electron-positron and hadron colliders is analysed. 
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1. Introduction 



A complete picture for both precise tests of the Standard Model [0] and a search of 
effects from new physics supposes a direct measurement of the three-boson electroweak 
vertex, searches of higgs particles 0], the supermultiplets etc. at colliders of super 
high energies (LEP200, LHC) as well as a study of the CP violation and a measurement 
of the fundamental parameters of the electroweak theory (first of all, in the heavy quark 
sector). 

In the nearest decade, the centre of efforts, directed to the realization of this program, 
will certainly be in the field of the heavy quark physics at both the running colliders (LEP 
and Fermilab) and the B meson factories, being planned in SLAC, KEK and at HERA-B. 
In this case, the extraction of effects, related with high values of the energy scale, will 
be essentially determined by an accuracy of the theoretical and empirical knowledge on 
mechanisms of the quark interactions at non high energy and, first of all, about effects, 
caused by the QCD dynamics @]. Therefore, the experimental researches of processes 
with the heavy c-, 6-, t-quarks take a special importance. 

The presence of the small parameter Aqcd/^^q, where Aqcd is the scale of the quark 
confinement and mg is the heavy quark mass, has allowed one to develop the powerful 
tools for the study of the QCD dynamics in the heavy quark interactions, such methods as 
the phenomenological potential models [p|-[T(]|], the QCD sum rules [^^-[T^ and Effective 
Heavy Quark Theory (EHQT) , that is successfully applied for the study of hadrons, 
containing a single heavy quark. Thus, the investigation of the processes with the heavy 
quarks allows one to extract and to study nonperturbative QCD effects, causing the 
quark hadronization, by means of the use of the heavy quark as the "marked" atoms. A 
successful realization of such program of studies becomes possible due to the progress in 
the experimental technique of the detecting and the identification of particles (mainly, it 
is related with the invention and the improvement of the vertex detectors, allowing one 
to observe the heavy quark particles due to its running gap from the primary vertex of 
the interaction). 

Among the heavy quarkonia {QQ'), the {be) system with the open charm and beauty 
takes a particular place. In contrast to the hidden charm (cc) and beauty {bb) families, 
studied in details experimentally [TH] and being quite accurately described theoretically 
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T^,[TB|,|T^, the heavy quarkonium (be), the family of mesons, has some specific pro- 



duction and decay mechanisms and the spectroscopy, whose study allows one to extend 
and to enforce the quantitative understanding of the QCD dynamics as well as to step 
forward in the study of the most important parameters of the electroweak theory. 

From the spectroscopy viewpoint, the (be) system is the heavy quarkonium, whose 
spectrum can be quite reliably calculated in the framework of the QCD-motivated non- 
relativistic potential models as well as in the QCD sum rules, (be) is the only system, 
composed of two heavy quarks, where the description of this system mass spectrum can 
test the selfconsistency for the potential models and the QCD sum rules, whose parameters 
(the quark masses, for instance) have been fixed from the fitting of the spectroscopic data 
on the charmonium and bottomonium. Thus, the study of the family spectroscopy can 
serve for the essential improvement of the quantitative characteristics of the quark models 
and the QCD sum rules, which are intensively applied in other fields of the heavy quark 
physics (for example, when one extracts values of elements in the matrix of mixings of 
the heavy quark weak charged currents and one estimates contributions, interfering with 
the effects of the CP invariance violation, in the heavy hadron decays ||18|| ). 

Moreover, there is a problem of the precise description of the P-wave level splittings 
in the charmonium and bottomonium, when the experimental measurement has found an 
essential deviation from the values, which have been expected in some well-acknowledged 



quark models ||19[. The study of the B^ meson family can help in a solution of this 
problem. 

In addition, the (be) system is interesting due to that it allows one, in a new way, 
to use the phenomenological information, obtained from the detailed experimental study 
of the charmonium and bottomonium. So, for example, (be) takes an intermediate place 
between the charmonium and bottomonium in respect to both the system level masses 
and the values of average distances between the heavy quarks. As has been clarified, in 
the region of the average distances in the (cc) and (bb) systems, the heavy quark potential 
possesses the simple scahng properties P, |2(]| , |57| , which state that the kinetic energy of 
the heavy quarks is practically a constant value, independent of the quark flavours and 
the excitation level in the heavy quarkonium system. Furthermore, this leads to that 
the heavy quarkonium level density (the distance between the nL and n'L levels) does 
not depend on the flavours of quarks, composing the heavy quarkonium. This regularity 
is quite accurately valid empirically for the (cc) and (bb) systems and it can be used in 
the framework of the QCD sum rules, where a scaling relation, connecting the leptonic 
constants of the S'-wave levels in the different quarkonia [21,2^, is derived. 

Further, having no strong and electromagnetic annihilation channels of decays, the 
excited [be) system levels, being below the threshold of the decay into the BD meson pair, 
will decay into the lightest basic pseudoscalar state 5^(0") due to the radiative cascade 
transitions into the underlying levels. Therefore, the widths of the electromagnetic (7) 
and hadronic (vrvr, 77, ...) radiative decays of the given excitation into the other levels will 
compose its total width. As a result, the total widths of the excited levels in the (be) 
system turn out two orders of magnitude less than the total widths of the charmonium and 
bottomonium excited levels, for which the annihilation channels are essential. Moreover, 
maybe, the data on the radiative hadronic decays in the {be) family give a possibility 
for one to solve some problems on the theory of the hadronic transitions in the heavy 
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quarkonia (for example, the problem on the anomalous distribution over the tttt pair 
invariant mass in the decay of T" — > Ttttt [p^ -p8|). 

Thus, on the one hand, the methods, applied in the heavy quark physics, are able quite 
reliably to point out the spectroscopic characteristics of the (6c) system for one to make 
a purposefully-directed experimental search of the given heavy quarkonium, and, on the 
other hand, the measurement of the spectroscopic data in the Be family would allow one 
to improve these methods approaches for the extraction of the fundamental parameters 
of the Standard Model from both the Be meson physics and the other fields of the heavy 
quark physics. 

Like the other mesons with the open flavour, the basic state of the Be meson family, 
the pseudoscalar meson i?+(0~), is the long-living particle, decaying due to the weak 
interaction and having the life time, comparable with the lifetimes of B and D mesons, so 
this feature essentially distinguishes Be from the heavy quarkonia rje and rji,. Therefore, the 
study of Be meson decays is the rich field of the heavy quark physics, where one extracts 
an important information about both the QCD dynamics and the weak interactions. The 
spectroscopic Be meson characteristics such as the leptonic constant, determining the 
width of the wave package of the (be) system in the basic state, essentially determine 
the description of the Be decay modes, in which some specific features and effects are 
observed. 

First of all, the presence of the valent heavy quark-spectator leads to a large probability 
for the Be decay modes with the heavy mesons in the final state, i.e. in the decays 
Be — > iljiVc) and Be — >■ B^J"^ p9|-p6|. The large ip particle yield is interesting, in addition, 
by that the ip particle has the perfect experimental signature in the leptonic decay mode. 

Further, in the consideration of the semileptonic B^ il){rie)l^v decays, the non- 
relativistic heavy quark motion inside the quarkonia leads to an essential effect, caused 
by large Coulomb-like corrections, which notably change the calculation results for these 
decays in the framework of the QCD sum rules |3^. Only the taking into the account 
these corrections makes the results of the QCD sum rules and the potential quark models 
to be consistent. 

Recently, considering the semileptonic transitions of the heavy quarks Q Q'lv in the 
framework of the Effective Heavy Quark Theory (EHQT) for hadrons with a single heavy 
quark (Qg, Qqq), one has stated the universal regularities [0], which serve, for example, 
for the model independent extraction of the Kobayashi-Maskawa matrix element value 
iVfecl- This universality in the limit of Kqcd/i^q — > is caused by the heavy quark 
flavour-independence of the light quark motion in the gluon field of the static source (the 
heavy quark), so that the wave functions of such hadrons are universal. In the case of 
the heavy quarkonium with two heavy quarks, the distances between the quarks depend 
on the values and ratios of its masses, i.e. the wave functions of the heavy quarkonia are 
not universal and depend on the quark flavours. However, in this case, one can neglect 
a low value of the spin-dependent splitting in the heavy quarkonium and suppose the 
wave functions of the nLj quarkonia to be J-independent. This fact finds the expression 
in an approximate spin-symmetry for the heavy quarks, so it puts some relations on the 



form factors of the weak semileptonic exclusive decays of Be Such relations for 

the form factors are universal and characteristic for the Be meson and reflect the high 
power of understanding the heavy quark decay dynamics, needing a direct experimental 
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verification. 

Considering the Be decays with the spectator 6-quark, one has particularly to note 
an essential role of the effects, caused by that the c-quark is not in the free state, but in 
the bound one. The decrease of the phase space for the c-quark decay within the heavy 
quarkonium makes the probability of the decay to be 40 % less than the probability in 
the D and Ds meson decays ||3^. The annihilation channel of the weak Be meson decay 



5^, allowing one to determine the value of the quark wave function at the origin |\E'(0)|^ 



acquires the important meaning, too. 

As in the case of the {be) system spectroscopy, the heavy quark theory is able to make 
the basic predictions on the mechanisms of the Be meson decays, whose characteristics 
measurement would allow one essentially to develop the methods of its description and 
also to use these methods for the precise investigations of the Standard Model as well as 
possible deviations from predictions of the latter. 

In the case of the Be meson production, a low value of the Kqcd/^q ratio and, 
hence, the low value of the quark-gluon coupling as ~ 1/ Inimq/ Kqcd) ^ 1 allow one 
to make the consideration of the pair production of the hh and cc quarks, from which the 
Be meson is formed, in the framework of the perturbative QCD theory, and also, in a 
way, to factorize contributions, caused by the perturbative production of heavy quarks 
and forthcoming nonperturbative binding of the latter into the heavy quarkonium. So, 
to calculate the cross sections of the S'-wave Be state production in the Z boson peak is 
enough to compute the matrix elements for the joint production of the hh and cc pairs 
in the colour-singlet state of the (he) pair with the fixed total spin of quarks {S = 0, 1), 
when the quarks, being bound into the meson, move with one and the same velocity, 
equal to the meson velocity. After that, one has to multiply these matrix elements by the 
nonperturbative factor, whose value is determined by the spectroscopic characteristics of 
the bound state (the quark masses and the leptonic constant, related with the probability 
of the observation of quarks at the zero distance between them in the bound state) [P^-p7| . 
The last notion is caused by that the characteristic virtualities of heavy quarks inside the 
heavy quarkonium are much less than its masses, since the heavy quarks inside the bound 
states are moving nonrelativistically, otherwise the quark virtualities in its production are 
of the order of its masses. Therefore, considering the Be production, one can assume, 
that, inside the meson, the h- and c-quarks are close to the mass shell and practically at 
rest in respect to each other. Thus, after the extraction of the nonperturbative factor, 
the analysis of the Be heavy quarkonium production is determined by the consideration 
of the matrix elements, calculated in the perturbation theory of QCD. 

Note first of all, that the necessity of the two pair production of heavy quarks in the 
electromagnetic and strong processes for the Be yield leads to that the leading order of 
the perturbative QCD has an additional factor of the suppression ~ a| in respect to the 
leading order of the perturbation theory for the production of the single flavour heavy 
quarks, for example, the hh pair (see Figures 0, ||), so a{Bc) / a{bb) ~ Q;||\E'(0)p/m|!. This 
causes the low yield of the Be mesons in respect to the B meson production. 

The analysis of the leading approximation in the perturbative QCD for the Be meson 
production allows one to derive a number of analytical expressions for the Be production 
cross sections [^,^ , where one has especially to stress the expressions for the functions of 
the heavy quark fragmentation into the heavy quarkonium in the scaling limit M^/ s — 0, 
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so these functions are determined by the values of a^, the quark masses and the leptonic 
constant of the meson [^^-Q. Thus, the fragment at ional Be production can be rehably 
described by the analytical expressions, so this opens new possibilities in the study of 
the QCD dynamics, essential in the complete picture of the heavy quark physics. As one 
can show, the fragmentational Be production certainly dominates in the Z boson decays 
so that it can be straightforwardly studied at the LEP facilities. Moreover, one can 
analytically study notable spin effects in the fragmentation into the vector B* meson , 
decaying electromagnetically B* Bc'j. 

In the hadronic Be production, patron processes at the energies, comparable with the 
Be mass, dominate, so that the processes, having the character of the fragmentational 
and also recombinational type P5|JiB[] (see Figure y), are essential. 

Further, the numerical estimates of the Be meson yield at the colliders LEP and 
Tevatron show that the fraction of the Be mesons in the production of the beauty hadrons 
is of the order of 10^^ p8|-^,^. This leads to that at the current experimental facilities, 
a quite large number of the Be mesons are being produced. 

Thus, one can point out the expected number of the Be mesons, being produced at 
different colliders, and the differential Be characteristics, whose experimental study would 
significantly clarify the picture of the QCD interactions of heavy quarks. 

A solution of the problem on the experimental discovery and study of the Be mesons 
is determined, first, by the theoretical description of the characteristics of the Be meson 
family (the spectroscopy, the production and decay mechanisms), so that the present 
review is devoted to this purpose. Second, this program is determined by the experimental 
methodics at the current detectors, so that the latter would allow one to observe the events 
with the Be production and decays, predicted by the theory. As for the second part of the 
problem, at present there is, as mentioned, the colossal progress, related with the use of 
the electronic vertex detectors, possessing the fast operation and allowing one to isolate 
the processes with the long living particles {B, Be, D) from the production processes 
(the technique of distinguishing the primary and secondary vertices) , and also accurately 
to reconstruct the decay vertices of the particles in space |Q. The presence of distinct 
signatures in the Be meson decays and the practical possibility for the registration of 
these decay modes have led to the real chance of the Be meson discovery at the LEP and 
Fermilab detectors as well as to the sharp rise of the theoretical interest to the {be) 
system. The latter has reflected in the achievement of a large number of the essential 
results in the consideration of the heavy quark interaction mechanisms at the example of 
Be mesons. So, the present paper is devoted to the review of these results. 



2. Spectroscopy of Be mesons 



Some preliminary estimates of the bound state masses of the (6c) system have been 
made in P,^, devoted to the description of the charmonium and bottomonium proper- 
ties, as well as in ref. |^2|. Recently in refs.||53| and the revised analysis of the Be 



spectroscopy has been performed in the framework of the potential approach and QCD 
sum rules. 

In the present section we consider the (be) spectroscopy with account of the change 
of the effective Coulomb interaction constant, defining spin- dependent splittings of the 
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quarkonium levels. We calculate the widths of radiative transitions between the levels and 
analyse the leptonic constant Jb^ in the framework of the QCD sum rules in the scheme, 
allowing one to derive scaling relation for the leptonic constants of the heavy quarkonia. 



2.1. Mass spectrum of Be mesons 

The Be meson is the heavy (be) quarkonium with the open charm and beauty. It 
occupies an intermediate place in the mass spectrum of the heavy quarkonia between the 
(cc) charmonium and the [bb) bottomonium. The approaches, applied to the charmonium 
and bottomonium study, can be used to describe the Be meson properties, and experi- 
mental observation of Be could serve as a test for these approaches and it could be used 
for the detailed quantitative study of the mechanisms of the heavy quark production, 
hadronization and decays. 

In the following we obtain the results on the Be meson spectroscopy. We will show 
that below the threshold for the hadronic decay of the (be) system into the BD meson 
pair, there are 16 narrow bound states, cascadely decaying into the lightest pseudoscalar 
B^{0~) state with the mass m(0~) ^ 6.25 GeV. 

2.1.1. Potential 

The mass spectra of the charmonium and the bottomonium are experimentally studied 
in details ]T3[ and they are properly described in the framework of phenomenological 
potential models of nonrelativistic heavy quarks P-^,p!0[] . To describe the mass spectrum 
of the (be) system, one would prefer to use the potentials, whose parameters do not depend 
on the flavours of the heavy quarks, composing a heavy quarkonium, i.e. one would use 
the potentials, which rather accurately describe the mass spectra of (cc) as well as (bb), 
with one and the same set of potential parameters. The use of such potentials allows 
one to avoid an interpolation of the potential parameters from the values, fixed by the 
experimental data on the (cc) and (bb) systems, to the values in the intermediate region 
of the (be) system. 



As it has been shown in ref. EO], with an accuracy up to an additive shift, the potentials. 



independent of heavy quark flavours [^|-^,10|, coincide with each other in the region of 



the average distances between heavy quarks in the (cc) and (bb) systems, so 

0.1 fm < r < 1 fm , (1) 

although those potentials have different asymptotic behaviour in the regions of very low 
(r — 0) and very large (r —>■ oo) distances. 

In Cornell model in accordance with asymptotic freedom in QCD, the potential has 
the Coulomb-like behaviour at low distances, and the term, confining the quarks, rises 
linearly at large distances 

Vc{r) = -^^ + -^+eo, (2) 
6 r 

so that 

as = 0.36 , 
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a = 2.34 GeV^^ , 
rric = 1.84 GeV , 

Co = -0.25 GeV. (3) 

The Richardson potential [0] and its modifications in refs. [|lOl and 0] also correspond 
to the behaviour, expected in the framework of QCD, so 

v«M . .f^,.^.i^ 1 



2^ 



(27r)3 3 llNc-2nf q'^\n{l + q^/A 

i!^ 1 ^ ( i ^ AV (4) 

[2ny 3 27 U'ln(l + gVA2) 27 ' ^' 



with 



A = 0.398 GeV. (5) 

In the region of the average distances between heavy quarks (|1]), the QCD-motivated 
potentials allow the approximations in the forms of the power (Martin) or logarithmic 
potentials. 

The Martin potential has the form ||^ 

Vuir) = -Cm + dM^^urf , (6) 

so that 



(7) 



Am 


= 1 GeV, 


k 


= 0.1, 


rrib 


= 5.174 GeV, 


rric 


= 1.8 GeV, 


Cm 


= 8.064 GeV, 




= 6.869 GeV. 



The logarithmic potential is equal to ^ 

VLir) = CL + dL In(Air) , (8) 

so that 

Al = 1 GeV, 
nib = 4.906 GeV, 

rric = 1.5 GeV, (9) 
cl = -0.6635 GeV, 
dL = 0.733 GeV. 

The approximations of the nonrelativistic potential of heavy quarks in the region of 
distances (|l|) in the form of the power (^) and logarithmic (H) laws, allow one to study its 
scaling properties. 
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In accordance with the virial theorem, the average kinetic energy of the quarks in the 
bound state is determined by the following expression 

1 rdV 

(T) ^ ^ . (10) 

Then, the logarithmic potential allows one to conclude, that for the quarkonium states 
one gets 

{Tl) = const (11) 
independently of the flavours of the heavy quarks, composing the heavy quarkonium, 

dL/2 = const ^ 0.367 GeV . 

In the Martin potential, the virial theorem ( p!0D allows one to obtain the expression 

(Tm) = ^ {cm + E) , (12) 

where E is the binding energy of the quarks in the heavy quarkonium. Phenomenolog- 
ically, one has \E\ <^ cm (for example, E{lS,cc) ~ —0.5 GeV), so that, neglecting the 
binding energy of the heavy quarks inside the heavy quarkonium, one can conclude that 
the average kinetic energy of the heavy quarks is a constant value, independent of the 
quark flavours and the number of the radial or orbital excitation. The accuracy of such 
approximation for (T) is about 10%, i.e. \AT/T\ ~ 30 ^ 40 MeV. 

From the Feynman-Hellmann theorem for the system with the reduced mass /z, one 

has 

f--^^, (13) 

djj, jj, 

and, in accordance with condition (pA]), it follows that the difference of the energies for 
the radial excitations of the heavy quarkonium levels does not depend on the reduced 
mass of the QQ' system 

E{n, fi) - E{n, fi) = E{n, fi') - E{n, fi') . (14) 

Thus, in the approximation of both the low value for the binding energy of quarks and 
the zero value for the spin-dependent splittings of the levels, the heavy quarkonium state 
density does not depend on the heavy quark flavours 

dn 

const . (15) 



dM, 



The given statement has been also derived in ref. ||21| using the Bohr-Sommerfeld quan- 
tization of the S-wave states for the heavy quarkonium system with Martin potential 

i- 

Relations (|1^-(|15]) are phenomenologically confirmed for the vector S-levels of the 
(bb), (cc), (ss) systems |]T^ (see Table 



Thus, the structure of the nonsplitted S-levels of the (be) system must repeat not only 
qualitatively, but quantitatively the structure of the S-levels for the (bb) and (cc) systems, 
with an accuracy up to the overall additive shift of masses. 
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Table 1. 



The mass difference (in MeV) for the two hghtest vector states of different heavy 
systems, AM = M{2S) - M{IS) 



System 


T 




Be 





AM 


563 


588 


585 


660 



Moreover, in the framework of the QCD sum rules, the universahty of the heavy quark 
nonrelativistic potential (the independence on the flavours and the scaling properties (pJ]), 
(0)' (HI)) allows one to obtain the scaling relation for the leptonic constants of the S-wave 
quarkonia [pl|] 

f 

— = const (16) 
independently of the heavy quark flavours in the regime, when 

\mQ-mQ,\ is restricted , 

i.e., when one can neglect the heavy quark mass difference. On the other hand, in the 
regime, when the mass difference is not low, one has 

S(4^)— 

where 



_ niQmQi 



rriQ + TTiQ' 

Consider the mass spectrum of the (be) system with the Martin potential [^. 

Solving the Schrodinger equation with potential (P) and the parameters (|^, one flnds 
the Be mass spectrum and the characteristics of the radial wave functions -R(O) and -R'(O), 
shown in Tables ^ and |^, respectively. 

The average kinetic energy of the levels, lying below the threshold for the (be) system 
decay into the BD pair, is presented in Table wherein one can see that the term, added 
to the radial potential due to the orbital rotation, 

m = ^, (18) 



Table 2. 

The energy levels of the (be) system, calculated without taking into account relativistic 
corrections, in GeV 



n 




52 




|55| 




54 




n 




521 


|55 






54 1 


n 


|52| 


|55 




|54| 


IS 


6.301 


6.315 


6.344 


2P 


6.728 


6.735 


6.763 


3D 


7.008 


7.145 


7.030 


2S 


6 


.893 


7.009 


6.910 


3P 


7.122 






7.160 


4D 


7.308 






7.365 


3S 


7.237 




7.024 


4P 


7.395 










5D 


7.532 
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Table 3. 

The characteristics of the radial wave functions Rnsi^) (in GeV^/^) and R'^p{0) (in 
GeV^/^), obtained from the Schrodinger equation 



n 


Martin 


H 




1.31 


1.28 


R2S{0) 


0.97 


0.99 


^2P(0) 


0.55 


0.45 




0.57 


0.51 



weakly influences the value of the average kinetic energy, and the binding energy for 
the levels with L 7^ is essentially determined by the orbital rotation energy, which is 
approximately independent of the quark flavours (see Table ^, so that the structure of the 
nonsplitted levels of the (be) system with L 7^ must quantitatively repeat the structure 
of the charmonium and bottomonium levels, too. 

2.1.2. Spin-dependent splitting of the (be) quarkonium 



In accordance with the results of refs. ||55|J56[] , one introduces the additional term to the 
potential to take into the account the spin-orbital and spin-spin interactions, causing the 
splitting of the nL levels (n is the principal quantum number, L is the orbital momentum), 
so it has the form 

4 1 L ■ S 4 2 

+ ^ as T- + ^ as ^ Sc ■ Sb 47r 5{r) (19) 

4 1 1 r 

+ «s (3(Sc ■ n) (Sfe • n) - Sc ■ Sfe) — , n = - . 

where V{r) is the phenomenological potential, conflning the quarks, the flrst term takes 
into account the relativistic corrections to the potential V{r)] the second, third and fourth 
terms are the relativistic corrections, coming from the account of the one gluon exchange 
between the b and c quarks; as is the effective constant of the quark-gluon interaction 
inside the (be) system. 

The value of the as parameter can be determined in the following way. 



Table 4. 

The average kinetic and orbital energies of the quark motion in the (be) system, in 
GeV 



nL 


IS 


2S 


2P 


3P 


3D 


(T) 


0.35 
0.00 


0.38 
0.00 


0.37 
0.22 


0.39 
0.14 


0.39 
0.29 
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The splitting of the S-wave heavy quarkonium {Q1Q2) is determined by the expression 

AM(n^) = as |i?„5(0)P , (20) 

9mim2 

where i?„5(0) is the value of the radial wave function of the quarkonium, at the origin. 
Using the experimental value of the S-state splitting in the (cc) system [IT3 

AM(15, cc) = 117 ± 2 MeV , (21) 

and the -Ris(O) value, calculated in the potential model for the (cc) system, one gets the 
model-dependent value of the q;s(V') constant for the effective Coulomb interaction of the 
heavy quarks (in the Martin potential, one has asiip) = 0.44). 

In ref.[Q the effective constant value, fixed in the described way, has been apphed to 
the description of not only the (cc) system, but also the (he) and (hb) quarkonia. 

In the present paper we take into account the variation of the effective Coulomb 
interaction constant versus the reduced mass of the system (/x) . 

In the one-loop approximation at the momentum scale p^, the "running" coupling 
constant in QCD is determined by the expression 



47r 



^sip') = ,,,,,,,,2 ^ ' (22) 



where b = 11 — 2nf/3, and nf = 3, when one takes into account the contribution by the 
virtual light quarks, < m^^. 

In the model with the Martin potential, for the kinetic energy of quarks (cc) inside ip, 
one has 

(Ti5(cc)) ^ 0.357 GeV , (23) 
so that, using the expression for the kinetic energy, 

{T) = f^, (24) 

one gets 

All 

"^(^^) = h^ fo/rr\ /A2 V ' (25) 

so that asiip) = 0.44 at 

Aqcd ~ 164 MeV . (26) 



Table 5. 



The average energy of the orbital motion in the heavy quarkonia, in the model with 
the Martin potential, in GeV 



System 


cc 


be 


bb 


AVi{2P) 


0.23 


0.22 


0.21 
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Table 6. 



The leptonic decay constants of the heavy quarkonia, the values, measured experi- 
mentally and obtained in the model with the Martin potential, in the model with the 
effective Coulomb interaction and from the scaling relation (SR), in MeV 



Model 


Exp. |T5| 


Martin 


Coulomb 


SR 


u 


410 ± 15 


547 ± 80 


426 ± 60 


410 ± 40 






510 ±80 


456 ± 70 


460 ± 60 


fr 


715 ±15 


660 ± 90 


772 ± 120 


715 ± 70 



As it has been noted in the previous section, the value of the kinetic energy of the 
quark motion weakly depends on the heavy quark flavours, and it, practically, is constant, 
and, hence, the change of the effective as coupling is basically determined by the variation 
of the reduced mass of the heavy quarkonium. In accordance with eqs.(25)-(p6D and Table 
^, for the (be) system one has 

TiL IS 2S 2P 3P 3D 
as 0.394 0.385 0.387 0.382 0.383. 

Note, the Martin potential leads to the -Ri5(0) values, which, with the accuracy up 
to 15 ^ 20%, agrees with the experimental values of the leptonic decay constants for the 
heavy (cc) and (bb) quarkonia. The leptonic constants are determined by the expression 



r{QQ ^ /+/- 

where eg is the heavy quark charge. 
In the nonrelativistic model one has 



— e' a' ^ 

3 ""Q "em ^ 



(27) 



QQ 



QQ 



\ 



TlM, 



RisiO) 



(2^ 



QQ 



For the effective Coulomb interaction of the heavy quarks in the basic IS-state one has 



M x3/2 

2 [- /. as) 



(29) 



One can see from Table that, taking into account the variation of the effective as- 
constant versus the reduced mass of the heavy quarkonium (see eq.(^)), the Coulomb 
wave functions give the values of the leptonic constants for the heavy IS-quarkonia, so 
that in the framework of the accuracy of the potential models, those values agree with the 
experimental values and the values, obtained by the solution of the Schrodinger equation 
with the given potential. 

The consideration of the variation of the effective Coulomb interaction constant be- 
comes especially essential for the T particles, for which ^^(T) ^ 0.33 instead of the fixed 
value as = 0.44. 
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Thus, calculating the splitting of the (be) levels, we take into account the as depen- 
dence on the reduced mass of the heavy quarkonium. 

As one can see from eq.(^), in contrast to the LS-coupling in the (cc) and (bb) systems, 
there is the jj-coupling in the heavy quarkonium, where the heavy quarks have different 
masses (here, LSc is diagonahzed at the given Jc momentum, (Jc = L + Sc, J = Jc + Sfe), 
J is the total spin of the system). We use the following spectroscopic notations for the 
splitted levels of the (be) system, - n^^^Lj. 

One can easily show, that independently of the total spin J projection one has 

l^^+^W) = |J = L + 1, ^ = 1), 

I'^'-'^L-i) = |J = L-1, 5 = 1) , (30) 



' ' \2L + r ' ^ V2L + 1' ' ' 



I'^-'Ll) = J^^\J = L, S = 1)-J^—\J = L, S = 0) , 

where | J, S) are the state vectors with the given values of the total quark spin S = Sc + S;,, 
so that the potential terms of the order of l/mjrii,, I /ml lead, generally speaking, to the 
mixing of the levels with the different Jc values at the given J values. The tensor forces 
(the last term in eq.(0)) are equal to zero at L = or S" = 0. 

To calculate values of the level shifts, appearing due to the spin-spin and spin-orbital 
interactions, one has to take the averaged expression (|T9[) over the wave functions of the 
corresponding states. 

The averaging over the angle variables can be performed in the following standard 
way. Let us represent the matrix element of the unit vector n = r/r pair in the form 

(L, mlnPn'^m m') = a(U'L'' + L''L^')w + M^'^^w , (31) 

where L are the orbital momentum matrices in the corresponding irreducible representa- 
tion. 

From the conditions of the normalization of the unit vector, {n'Pn'^)6^'^ = 1, the or- 
thogonality of the radius-vector to the orbital momentum, n^L^ = 0, the commutation 
relations for the angle momentum, [L^; L"^] = ie^'^^'Li, one finds the values of constants a 
and b in eq. (|3l|) 

1 



4L2 -3 
4L2 -3 ■ 



(32) 
(33) 



Note further, that from the condition for the quark spins SqSq + SqSq = it follows, 
that 

3(n%'? - -6P'')SPS^^ = -(n^n^ - -6P'')SpS'^ . (34) 



Thus, (see also ref . ||57| ) 



(6(n%^ - Is^S^^S',) = -^J^(6(LS)2 + 3(LS) - 2VS') . (35) 
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Fig. 1. The mass spectrum of the (be) system with account of sphttings. 



Using eqs.(^Dp, (|53|), for the level shifts, calculated in the perturbation theory at S' = 1 
one gets the following formulae 

2 



AEn3p2 
AEniPf, 



-as 



as 



Smcrrib 
2 



6 



\Rnsm' 



as 



n^D3 



-as 
as 



brricmh r-^ 4 m^J rdr 

4 , 1 . 1/1 1 \ ^ dV{r 

b 



AE, 



52 
21mcmb 
92 



-as 



^r3^ 2 {ml ^ ml) ^ 
r'^ 2 \m^ mi J 



rdr 
dy(r; 



3 

4 vrri 



rdr 3 
<W{r 



3 r'^ 
8 1 , 



rdr 



(36) 
(37) 
(38) 
(39) 
(40) 
(41) 



21m 

where Rns{0) are the radial wave functions at L = 0, (...) denote the average values, 
calculated under the wave functions RnL{f)- The mixing matrix elements have the forms 



( ^Pi\AE\ -'Pi) 



-as 



9mcmb 



(4) + 
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Table 7. 



The masses (in GeV) of the lightest pseudoscalar Be and vector B* states in different 
models (* is the present paper) 



State 


* 




N 




m 


11 


i 


m 


mm 


0- 


6.253 


6.249 


6.314 


6.293 


6.270 


6.243 


6.246 




6.317 


6.339 


6.354 


6.346 


6.340 


6.320 


6.319 


State 


m 


m 




IP 


H 






in 




m 




0- 


6.264 


6.320 


6.256 


6.276 


6.286 




6.255 


1- 


6.337 


6.370 


6.329 


6.365 


6.328 


6.320 


6.330 



^Pi\AE\ ^Pi) 



( ^Pi\AE\ ^Pi) 



/ 1 



\4m^ 



-as 



12ml 



)<■ 



dV(r) 8 1 , 
rdr 3 r-^ 



4 



1 



1 



Til) ^' 



2ml 

2v^ 1 

V2 dVjr) ^ 8 
rdr 3 



dl^(r) 8 1 , 
3 r-^ 



rdr 



6rri^ 



1. 



4 



-as 



2m^ 

-"5 



1 1 \ dV{r) 

5ml J rdr 



1 



3 9 



Am^. 
as 



+ 



20m? 



dF(r) 8 



rdr 



2v^ /Iv 



Ibmcrrib 
V6 , dV{r 



lOmg 



1. 

rdr 3 r'^ 



+ 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



As one can see from eq.(^), the S-level splitting is essentially determined by the 
|-Rn5(0)| value, which can be related to the leptonic decay constants of the S-states (0~, 
1^). Section 2.3 is devoted to the calculation of these constants in different ways. We 
only note here, that with enough accuracy, the predictions of different potential models 
on the |-Ri5(0)| value are in agreement with each other as well as with predictions in other 
approaches. 

For the 2P, 3P and "iD levels, the mixing matrices of the states with the total quark 



15 



spin S = 1 and S = have the forms 

|2P, 1'+) = 0.294|^= 1) +0.956|^ = 0) , (48) 
|2P, 1+) = 0.956|5 = 1) -0.294|5 = 0) , (49) 

so that in the 1"*" state the probabihty of the total quark spin value S* = 1 is equal to 

Wi{2P) = 0.913 . (50) 



For the 3P level one has 



|3P, 1'+) = 0.371|S = 1) +0.929|S = 0) , (51) 
|3P, 1+) = 0.929|5 = 1) -0.371|5 = 0) , (52) 



so that 



wi{3P) = 0.863 , (53) 

For the 3-D level one gets 

|3P>, 2'-) = -0.566|^ = 1) + 0.825|5 = 0) , (54) 

|3P>, 2") = 0.825|S = 1) +0.566|5 = 0) , (55) 

so that 

W2i3D) = 0.680 . (56) 

With account of the calculated splittings, the Be mass spectrum is shown in Figure ^ 
and Table 

The masses of the mesons have been also calculated in papers of ref. [^. 

As one can see from Tables ^ and ^ the place of the IS-level in the {be) system 
(m(15') ~ 6.3 GeV) is predicted by the potential models with the rather high accuracy 
Am(lS') ~ 30 MeV, and the IS-level splitting into the vector and pseudoscalar states is 
about m(l") - m(0") ^ 70 MeV. 

2.1.3. Be meson masses from QCD sum rules 

Potential model estimates for the masses of the lightest (be) states are in agreement 
with the results of the calculations for the vector and pseudoscalar (be) states in the 
framework of the QCD sum rules [^,^,0, where the calculation accuracy is lower, 



than the accuracy of the potential models, because the results essentially depend on both 
the modelling of the nonresonant hadronic part of the current correlator (the continuum 
threshold) and the parameter of the sum rule scheme (the moment number for the spectral 
density of the current correlator or the Borel transformation parameter). 



m*^(0~) ^ m^^(r) ^ 6.3 ^6.5 GeV . (57) 

As it has been shown in [^, for the lightest vector quarkonium, the following QCD 
sum rules take place 

IVIy — Q TT J Si S — Q 
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Table 8. 

The masses (in GeV) of the bound (be) states below the threshold of the decay into 
the BD meson pair (* is the present paper) 



State 


* 


M 




0] 




I'So 


6.253 


6.264 


6.314 


l^Si 


6.317 


6.337 


6.355 


2^So 


6.867 


6.856 


6.88 


9 


2^Si 


6.902 


6.899 


6.917 


2'Po 


6.683 


6.700 


6.72 


8 


2P 1+ 


6.717 


6.730 


6.760 


2P 1'+ 


6.729 


6.736 






2 P2 


6.743 


6.747 


6.773 


S^Po 


7.088 


7.108 


7.134 


3P 1+ 


7.113 


7.135 


7.159 


3P 1'+ 


7.124 


7.142 








7.134 


7.153 


7.166 


3D 2- 


7.001 


7.009 








7.007 


7.005 






33^1 


7.008 


7.012 






3^ 2'- 


7.016 


7.012 







where fv is the leptonic constant of the vector (be) state with the mass My, 

tfvMve^^expitpx) = {0\Ux)\V{p, X)) , (59) 
J^(x) = c(x)7^6(x) , (60) 

where A, p are the B* polarization and momentum, respectively, and 

J d'x exp {zqx){0\TJ,ix)MO)\0) = i^-g,, + H^^^ + q,q, nf ° , (61) 

nQ^°(g2) = nQC°(P^'^*) + ng^°(°°°P^'^*)(g2) , (62) 

j^QCD(nonpert)^^2) ^ Y.Q{q^)0\ (63) 

where are the vacuum expectation values of the composite operators such as {rmpip), 
{as G^j,), etc. The Wilson coefficients are calculable in the perturbation theory of QCD. 
Si = {rric + rrif,)^ is the kinematical threshold of the perturbative contribution. My > Sj, 
Sth is the threshold of the nonresonant hadronic contribution, which is considered to be 
equal to the perturbative contribution at s > Sth- 

Considering the respective correlators, one can write down the sum rules, analogous 
to eg. (|58|) , for the scalar and pseudoscalar states. 

One believes that the sum rule ( ^8]) must rather accurately be valid at q^ < 0. 

For the n-th derivative of eq. (0) at = one gets 

^ i *I^J»£) , tf -^nr"-%^^ . (64) 
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so, considering the ratio of the n-th derivative to the n + 1-th one, one can obtain the 
value of the vector B* meson mass. The calculation result depends on the n number in 
the sum rules (|64D, because of taking into the account both the finite number of terms in 
the perturbation theory expansion and the restricted set of composite operators. 

The analogous procedure can be performed in the sum rule scheme with the Borel 
transform, leading to the dependence of the results on the transformation parameter. 

As one can see from eg. ([6^) , the result, obtained in the framework of the QCD sum 
rules, depends on the choice of the values for the hadronic continuum threshold energy and 
the current masses of quarks. Then, this dependence causes large errors in the estimates 
of the masses for the lightest pseudoscalar, vector and scalar (be) states. 

Thus, the QCD sum rules give the estimates of the quark binding energy in the 
quarkonium, and the estimates are in agreement with the results of the potential models, 
but sum rules involve a considerable parametric uncertainty. 

2.2. Radiative transitions in the Be family 

The Bf. mesons have no annihilation channels for the decays due to QCD and electro- 
magnetic interactions. Therefore, the mesons, lying below the threshold for the B and 
D mesons production, will, in a cascade way, decay into the 0~(15') state by emission of 
7 quanta and tt mesons. Theoretical estimates of the transitions between the levels with 
the emission of the vr mesons have uncertainties, and the electromagnetic transitions are 
quite accurately calculable. 

2.2.1. Electromagnetic transitions 



The formulae for the radiative El-transitions have form 17,68 



T{nPj^n^S^ + ^) = ^a,^Qlf,u' l\nP;nS)wj{nP) , 

T{nPj^n'So + l) = ^ a,^ Qls l\nP;nS) (1 - WjinP)) , 

Tin's, ^nPj + ^) = ^a^^QlsUj' l\nS;nP){2J + l)wj{nP) , (65) 

T{n'So^nPj + ^) = ^ a,^ Ql^ PinS;nP) {2J + 1) {1 - wjinP)) , 

y 

T{nPj riDj, + 7) = 1 a,^ Ql^ l\nD; nP) (2J' + 1) wj{nP))wj:{nD)Sjj. , 
T{nDj nPj. + 7) = A q2^ ^3 j2^^^. -p) ^^j, ^ ^) wr{nP))wj{nD)Srj , 

where uj is the photon energy, aem is the electromagnetic fine structure constant. 
In eq.(|65D one uses 

m^Ql - rribQc 

Qcs = ; , (06) 

rric + nih 

where Qcfi are the electric charges of the quarks. For the B^ meson with the parameters 
from the Martin potential, one gets Qeff = 0.41. 
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wjinL) is the probabihty that the spin S" = 1 in the nL state, so that Wo{nP) = 
W2{nP) = 1, WiijiD) = w^ijiD) = 1, and the wi{nP), W2{nD) values have been presented 
in the previous section (see eqs.(pO|), (p3[), (0) ). 

The statistical factor Sjji takes values [B3] 



J 


J' 


Sjj' 





1 


2 


1 


1 


1/2 


1 


2 


9/10 


2 


1 


1/50 


2 


2 


9/50 


2 


3 


18/25 



The I{nL] nL') value is expressed through the radial wave functions, 

I{nL-nL') = \ j RnL{r)RnL'{r)r^^r\ . (67) 

For the set of the transitions one obtains 

/(15, 2P) = 1.568 GeV"^ , /(15, 3P) = 0.255 GeV"^ , 

I{2S,2P) = 2.019 GeV-\ /(25, 3P) = 2.704 GeV"\ (68) 

I{3D, 2P) = 2.536 GeV"^ , I{3D, 3P) = 2.416 GeV"^ . 

For the dipole magnetic transitions one has |]5|, p!7| , |68[| 

16 

nn'S, ^ n'Sf + 7) = y f^ls (2/ + 1) , (69) 

where 

Aif = J Rnsir)Rnsir)3oiujr/2y dr , 

and 

Ateflf = ^ {Qcrrib - Qimc) . (70) 

2 2mcmh 

Note, in contrast to the if) and T particles, the total width of the B* meson is equal to 
the width of its radiative decay into the Bc{0~) state. 

The electromagnetic widths, calculated with accordance of eqs.(^),(^), and the fre- 
quencies of the emitted photons are presented in Tables ^ |l^, |Tl|. 

Note, EO-transitions with the conversion of virtual 7-quantum into the lepton pair can 
take place. Moreover, due to the tensor forces, the states with J > and S = 1 can, 
in addition to the L-wave, have the admixture of |L ± 2|-waves, giving the quadrupole 
moment to the corresponding states and causing the E2-transitions. However, the men- 
tioned transitions are suppressed by the additional factor aem in the first case, and by the 
small value of amplitude, determining, say, the probability of the admixture appearance 
of the D-wave in the l~{nS) state. 

Thus, the registration of the cascade electromagnetic transitions in the (be) family can 
be used for the observation of the higher (be) excitations, having no annihilation channels 
of the decays. 
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The energies (in MeV) and widths (in keV) of the electromagnetic El-transitions in 
the (be) family (* is the present paper) 



Transition 


UJ 


r 


* 




rg 


2P2 - 


►1^1 + 7 


426 


102.9 


112.6 


2Po - 


► 15i + 7 


366 


65.2 




79.2 


2P 1'+ 


^1^1+7 


412 




.1 




0.1 


2P 1+ 


^ ISi + 7 


400 


77.5 




99.5 


2P 1'^ 


^ 15o + 7 


476 


131 


.1 


56.4 


2P 1+ 


^1^0 + 7 


464 


11.6 


0.0 


3P2 - 


► ISi + 7 


817 


19.2 




25.8 


3Po - 


► ISi + 7 


771 


16.1 




21.9 


3P 1'^ 


^1^1+7 


807 




.5 


2.1 


3P 1+ 


^ ISi + 7 


796 


15.S 




22.1 


3P 1'^ 


^1^0 + 7 


871 


20.1 






3P 1+ 


^1^0 + 7 


860 




.1 






3P2 - 


> 25i + 7 


232 


49.4 


73.8 


3Po - 


> 25i + 7 


186 


25.^ 




41.2 


T~i 1 f-\- 

3P 


^ 2Si + 7 


222 


c 


.9 


5.4 


3P 1+ 


^25i + 7 


211 


32.1 




54. d 


3P 1'+ 


^25o + 7 


257 


55 


.0 




3P 1+ 


^25o + 7 


246 




.1 




2^1 ^ 


2P2 + 7 


159 


14.5 




17.7 


2^1 ^ 


2P0 + 7 


219 


7.7 


7.8 


2^1 


2P 1'+ + 7 


173 


1 


.0 


0.0 


2Si ^ 


2P 1+ + 7 


185 


12.5 




14.5 




2P 1'+ + 7 


138 


15.9 


5.2 


2So^ 


2P 1+ + 7 


150 


1 


.9 


0.0 



2.2.2. Hadronic transitions 

In the framework of QCD the consideration of the hadronic transitions between the 
states of the heavy quarkonium family is built on the basis of the multipole expansion for 
the gluon emission by the heavy nonrelativistic quarks [Q, with forthcoming hadroniza- 
tion of gluons, independently of the heavy quark motion. 

In the leading approximation over the velocity of the heavy quark motion, the action, 
corresponding to the heavy quark coupling to the external gluon field, 

S-.nt = -9 j A'x Al{x) ■ . (71) 

can be expressed in the form 

5int = 9 j^t r^El{t,^) ^^^{v)^\v) K{sn,f) d'r , (72) 
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The energies (in MeV) and widths (in keV) of the electromagnetic El-transitions in 
the (6c) family (* is the present paper) 



Transition 


UJ 


r 


* 






3P2 ^ 


3^1 + 7 


126 


0.1 




0.2 


3P2 ^ 


3D 2'- + 7 


118 


0.5 






3P2 ^ 


3P' 2" -h 7 


133 


1 






3.2 


3P2 ^ 


bp's + 7 


127 


10.9 


17.8 


3Po ^ 


3^1 + 7 


80 


a 


.2 




6.9 


3P r+ 


^ 3Pi + 7 


116 


0.3 




0.4 


3P 1+ 


^ 3Pi + 7 


105 


1 


.6 


0.3 


3P 1'+ 


3D 2'- + 7 


108 


a 


.5 




" 


3P 1+ 


^ 3P 2- +7 


112 


a 


.9 


9.8 


3P 1'+ 


^ 3P 2" + 7 


123 


2 






11.5 


3P 1+ 


3D 2'- + 7 


97 


1 


.2 






3P)3 ^ 


2P2 + 7 


264 


76.9 


98.7 


3P'i 


2Po + 7 


325 


79.7 


88.6 


3Di 


2P 1'+ + 7 


279 


a 


.a 




0.0 


3Di 


2P I+ + 7 


291 


39.2 




49.3 


3Di 


2P2 + 7 


265 


2 


.2 




2.7 


3D 2'- 


^ 2P2 + 7 


273 


6.8 






3D 2'- 


^ 2P2 + 7 


258 


12.2 


24.7 


3D 2'- 


2P l'+ + 7 


287 


46.0 


92.5 


3D 2'- 


^ 2P I+ + 7 


301 


25.0 




3D 2- 


2P 1'+ + 7 


272 


15 


.4 


0.1 


3D 2- 


^ 2P 1+ + 7 


284 


44.6 


88.8 



where \l/n(i') is the wave function of the quarkonium, emitting gluon, (r) is the wave 
function of the colour-octet state of the quarkonium, K{sn, f) corresponds to the spin fac- 
tor (in the leading approximation, the heavy quark spin is decoupled from the interaction 
with the gluons). 

Then the matrix element for the El-El transition of the quarkonium nLj — > n'L'j, -\-gg 
can be written in the form 

M{nLj ^ n'L'j, + gg) = A^as E^E^ ■ 

I dVd\' Gt,s„ir,r') vl/„^^(r)vi/„,i, , (r'), (73) 

where (j^^^^(r, r') corresponds to the propagator of the colour-octet state of the heavy 
quarkonium 

where H^q is the hamiltonian of the coloured state. 
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The energies (in MeV) and widths (in keV) of the electromagnetic Ml-transitions in 
the (6c) family (* is the present paper) 



Transition 


UJ 




r 


|53 


2Si ^ISo + i 
2S'o^lSi+7 
15*1 ISq + 7 
2^1 ^ 2So + 7 


649 
550 
64 
35 


0.098 
0.096 
0.060 
0.010 


0.123 
0.093 
0.135 
0.029 



One can see from eq.(73), that the determination of the transition matrix element 
depends on both the wave function of the quarkonium and the hamiltonian Hqq. Thus, 
the theoretical consideration of the hadronic transitions in the quarkonium family is model 
dependent. 

In a number of papers of ref. for the calculation of the values such as ([73[), the 
potential approach has been developed. 

In papers of ref. |^ it is shown that nonperturbative conversion of the gluons into 
the TT meson pair allows one to give a consideration in the framework of the low-energy 
theorems in QCD, so that this consideration agrees with the papers, performed in the 
framework of PC AC and soft pion technique ||2^ . 

However, as it follows from eq. (|73|) and the Wigner-Eckart theorem, the differential 



width for the El-El transition allows the representation in the form |24 



dr 

dm^ 



[riLj n'L'j, + h) = (2 J' 




AkiL,L') 



(75) 



where is the invariant mass of the light hadron system h, {} are 6j-symbols, Ak{L, L') 
is the contribution by the irreducible tensor of the rang, equal to /c < 3, s is the total 
quark spin inside the quarkonium. 

In the limit of soft pions, one has Ai{L, V) = 0. 

From eqs.([73|), ( [75| ) it follows, that, with the accuracy up to the difference in the 
phase spaces, the widths of the hadronic transitions in the {QQ) and {QQ') quarkonia 
are related to the following expression 



(76) 



r(QQ) {r^iQQ)y 



Then the experimental data on the transitions of tp' — > J/ip+Tni, T' — > T+tttt, iIj{3770) ^ 
J/ip + uTT ^\ allow one to extract the values of Ak{L, L') for the transitions 25" — >■ IS' + tttt 
and 3^ ^ 1^ + TTTT [0 



The invariant mass spectrum of the tt meson pair has the universal form |25,26[ 

1 dr „ Ik. 



B 



r dm M2 
where x = m/2m^, Ik^^-jrl is the tttt pair momentum. 



(2x2 - 1) - 1 



(77) 
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The estimates for the widths of the hadronic transitions in the (be) family have been 



made in ref.[£3|. The hadronic transition widths, having the values comparable with the 
electromagnetic transition width values, are presented in Table |T^. 

The transitions in the (6c) family with the emission of f] mesons are suppressed by the 
low value of the phase space. 

Thus, the registration of the hadronic transitions in the (be) family with the emission 
of the TT meson pairs can be used to observe the higher 2S- and 3D-excitations of the basic 
state. 

2.3. Leptonic constant of Be meson 

As we have seen in Section 2.1, the value of the leptonic constant of the meson 
determines the splitting of the basic IS-state of the (be) system. Moreover, the higher 
excitations in the (be) system transform, in a cascade way, into the lightest 0~ state of Be, 
whose widths of the decays are essentially determined by the value of Jb^, too. In the quark 
models |]69|-[7T1, used to calculate the weak decay widths of mesons, the leptonic constant. 



as the parameter, determines the quark wave package inside the meson (generally, the 
wave function is chosen in the oscillator form), therefore, the practical problem for the 
extraction of the value for the weak charged current mixing matrix element |Vf,c| from the 
data on the weak B^ decays can be only solved at the known value of Jb^- 

Thus, the leptonic constant Jb^ is the most important quantity, characterizing the 
bound state of the (be) system. 

In the present Section we calculate the value of /^^ in different ways. 

To describe the bound states of the quarks, the use of the nonperturbative approaches 
is required. The bound states of the heavy quarks allow one to consider simplifications, 
connected to both large values of the quark masses AqcD/fng <^ 1 and the nonrelativistic 
quark motion v —>■ 0. Therefore the value of Jb^ can be quite reliably determined in the 
framework of the potential models and the QCD sum rules |ll| . 



Table 12. The widths (in keV) of the radiative hadronic transitions in the (be) family 



Transition 


r m 


25*0 15*0 + TTTT 


50 


2^1 ^ 15i + irn 


50 


3^1 ^ ISi + TTTT 


31 


ISi + TTTT 


32 


3^3 -> 1^1 + vrvr 


31 


3D2 ISo + TTTT 


32 
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Table 13. 



The leptonic Be meson constant (in MeV), calculated in the different potential mod- 
els (the accuracy ~ 15%) 



Model 


Martin 


Coulomb 


i 


m 


m 




m 




510 


460 


570 


495 


410 


600 


500 



2.3.1. Jb^ from potential models 



In the framework of the nonrelativistic potential models, the leptonic constants of the 
pseudoscalar and vector mesons (see eqs.(^9D, (|60D ) 



{0\cixh,bix)\B:ip,e)) 
{0\c{xh,^^b{x)\B,{p)) 

are determined by expression (p8[) 



ifv My exp{ipx) , 
ifp exp{ipx) , 



fv = fp 



Bc{lS) 



-Risio) , 



(78) 
(79) 



(80) 



where -Ri5'(0) is the radial wave function of the IS' state of the (be) system, at the origin. 
The wave function is calculated by solving the Schrodinger equation with different po- 
tentials [p|-pt p!0| , ^] , in the quasipotential approach |7^ or by solving the Bethe-Salpeter 
equation with instant potential and in the expansion up to the second order over the 



quark motion velocity v/c |73,74 



The values of the leptonic meson constant, calculated in different potential models 
and effective Coulomb potential with the "running" as constant, determined in Section 
2.1, are presented in Table |13|. 

Thus, in the approach accuracy, the potential quark models give the Jb^ values, which 
are in a good agreement with each other, so that 



/P°* = 500 ±80 MeV. 



(81) 



2.3.2. /b^ from QCD sum rules 



In the framework of the QCD sum rules [0, expressions (|58D -(|64D have been derived 
for the vector states. The expressions are considered at < in the schemes of the 
spectral density moments (|6^ or with the application of the Borel transform [|ll| . As one 
can see from eqs.(|5^ - (p^, the result of the QCD sum rule calculations is determined 
not only by physical parameters such as the quark and meson masses, but also by the 
unphysical parameters of the sum rule scheme such as the number of the spectral density 
moment or the Borel transformation parameter. In the QCD sum rules, this unphysical 
dependence of the /^^ value is due to the consideration being performed with the finite 
number of terms in the expansion of the QCD perturbation theory for the Wilson coeffi- 
cients of the unit and composite operators. In the calculations, the set of the composite 
operators is also restricted. 
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Table 14. 



The leptonic Be constant (in MeV), calculated in the QCD sum rules (SR - the 
scaling relation) 



Model 


m 




m 




m 


in 


m 


m 


SRpT 






375 


400 


360 


300 


160 


300 


450 


460 



Thus, the ambiguity in the choice of the hadronic continuum threshold and the pa- 
rameter of the sum rule scheme essentially reduces the reliability of the QCD sum rule 
predictions for the leptonic constants of the vector and pseudoscalar states. 

Moreover, the nonrelativistic quark motion inside the heavy quarkonium v leads 
to the as/f -corrections, becoming the most important, to the perturbative part of the 
quark current correlators, where as is the effective Coulomb coupling constant in the 
heavy quarkonium. As it is noted in refs. [|Tli^[76l , the Coulomb Q;5/f-corrections can 



be summed up and represented in the form of the factor, corresponding to the Coulomb 
wave function of the heavy quarks, so that 

f(„) = l|^|l-exp(-W3«)|-'. (82) 

where 2v is the relative velocity of the heavy quarks inside the quarkonium. The expansion 
of the factor (^) in the first order over as/v 

F(v) « 1 + ^ . (83) 

gives the expression, obtained in the first order of the QCD perturbation theory [pT[] . 

Note, the as parameter in eq.(^2t) should be at the scale of the characteristic quark vir- 
tualities in the quarkonium (see Section 2.1), but at the scale of the quark or quarkonium 
masses, as sometimes one makes it thereby decreasing the value of factor (^2|). 

The choice of the as parameter essentially determines the spread of the sum rule 
predictions for the Jb^ value (see Table |l^) 

/If = 160 ^ 570 MeV . (84) 

As one can see from eq. (^4]) , the ambiguity in the choice of the QCD sum rule pa- 
rameters leads to the essential deviations in the results from the Jb^ estimates (|8l|) in the 
potential models. 

However, as it has been noted in Section 2.1, 

1) the large value of the heavy quark masses AgcD/fng ^ Ij 

2) the nonrelativistic heavy quark motion inside the heavy quarkonium u — > 0, and 

3) the universal scaling properties of the potential in the heavy quarkonium, when 
the kinetic energy of the quarks and the quarkonium state density do not depend on the 
heavy quark flavours (see eqs. ([lOD - (H)), 

allow one to state the scaling relation (^) for the leptonic constants of the S-wave 
quarkonia 



TT l~ = const 
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Indeed, at AqcD/f^Q ^ 1 one can neglect the quark-gluon condensate contribution, 
having the order of magnitude 0{l/mi,mc) (the contribution into the ip and T leptonic 
constants is less than 15%). At v —>■ one has to take into account the Coulomb-like 
a^/f -corrections in the form of factor (^2]), so that the imaginary part of the correlators 
for the vector and axial quark currents has the form 

QmUviq^) ^ 53mnp(g2) = ^ f^Y , (85) 



2 " V4/iy 
where 

- 1 - 4mbme V ^0 

— {nib — TTLcy ' 

Moreover, condition ([I5D can be used in the specific QCD sum rule scheme, so that this 
scheme excludes the dependence of the results on the parameters such as the number of 
the spectral density moment or the Borel parameter. 

Indeed, for example, the resonance contribution into the hadronic part of the vector 
current correlator, having the form 



n(res) {q^)= f^^J2 fvnMLSis ~ M^yJ , (86) 

J b Q ^ 



can be rewritten as 



nr'fe^)-/j4V^/?.<..^^i:«("-^-). (87) 

where n{s) is the number of the vector S-state versus the mass, so that 

n{ml) = k . (88) 
Taking the average value for the derivative of the step-like function, one gets 

nr'tf)H^i:«(B-*))/j^./?„,,^^. (89) 

and, supposing 

i^Eein-k))^!, (90) 

one can, in average, write down 

^m{U^^^'^'\q^)) = 53mn(Q^°)(g2) , (91) 

so, taking into the account the Coulomb factor and neglecting power corrections over 
l/niQ, at the physical points s„ = one obtains 

Mn U/i/ TT dn ' ^ ' 
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where one has supposed that 



mfe + me ^ Mb, , (93) 

fvn ~ fpn = fn ■ (94) 



Further, as it has been shown in Section 2.1, in the heavy quarkonium the value of 
dn/dM„ does not depend on the quark masses (see eq.(p!5D), and, with the accuracy up 
to the logarithmic corrections, as is the constant value (the last fact is also manifested 
in the flavour independence of the Coulomb part of the potential in the Cornell model). 
Therefore, one can draw the conclusion that, in the leading approximation, the right 
hand side of eq. (p2|) is the constant value, and there is the scaling relation ([T7|) pT[| . This 
relation is valid in the resonant region, where one can neglect the contribution by the 
hadronic continuum. 

Note, scaling relation (^) is in a good agreement with the experimental data on 
the leptonic decay constants of the and T particles (see Table 0), for which one has 



4/i/M = 1 IT 



The value of the constant in the right hand side of eq.(|T^ is in agreement with the 
estimate, when we suppose 

dM-T 1 

(— ^) ^ -[(Mt'-Mt) + (Mx"-MtO] , (95) 
dn 2 

and as = 0.36, as it is in the Cornell model. 

Further, in the limit case of B and D mesons, when the heavy quark mass is much 
greater than the light quark mass rriQ ^ rrig, one has 

and 

f-M=^-^^f.^. (96) 
TT dn 

Then it is evident that at one and the same fi one gets 

f M = const . (97) 



Scaling law (0) is very well known in EHQT |[T^ for mesons with a single heavy quark 



(Qq), and it follows, for example, from the identity of the B and D meson wave functions 
within the limit, when an infinitely heavy quark can be considered as a static source of 
gluon field (then eq.(^) follows from eq. (pO|) ). 

In our derivation of eqs.(|96D and (|97|) we have neglected power corrections over the 
inverse heavy quark mass. Moreover, we have used the notion about the light constituent 
quark with the mass equal to 

ruq ^ 330 MeV , (98) 

so that this quark has to be considered as nonrelativistic one v ^ 0, and the following 
conditions take place 

rriQ + mg ^ M^*^-) , rrig <^ niQ , (99) 

fv ~ fp = f. (100) 
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In agreement with eqs.(0) and (98), one finds the estimates^ 

fsM = 120 ± 20 MeV , (101) 
= 220 ± 30 MeV , (102) 

that is in an agreement with the estimates in the other schemes of the QCD sum rules 



11,12 



Thus, in the limits of 4/i/M = 1 and /i/M <^ 1, scahng relation ( |T7| ) is consistent. 
The Jb^ estimate from eq.(0) contains the uncertainty, connected to the choice of the 
ratio for the b- and c-quark masses, so that (see Table 

/b, = 460 ± 60 MeV . (103) 

In ref . the sum rule scheme with the double Borel transform was used. So, it allows one 
to study effects, related to the power corrections from the gluon condensate, corrections 
due to nonzero quark velocity and nonzero binding energy of the quarks in the quarkonium. 
Indeed, for the set of narrow pseudoscalar states, one has the sum rules 

E , ^%fu2 = - / ^ ^^np(.) + Ca{q') {^G^} , (104) 



where 



2, 1 g2 /3(3f2 + 1)(1 -t.2)2 v + i 9^;4 + 4^2^3 



and 



) = 192^;^ f [ 2^^ ^ V- ) ' 



Acting by the Borel operator LT-{—q^) on eq. ( |104| ), one gets 
where 



E , ^M-Mlr) = - Us ^mllp{s) exp(-.r) + C^(r) (^G^) , (107) 



nl \ ax J X 



n,x— >oo 



C'air) = U{~q') Cciq') . (109) 

For the exponential set in the left hand side of eq.( |107| ), one uses the Euler-MacLaurin 
formula 

^ (J'itv ^"P(-^^) = rdM,-^M,V|,exp(-M^) + 

n-l 

Y^Mtflu exp(-M2r) + ■ ■ • . (110) 



fc=0 



^In ref.|£T| the dependence of the S-wave state density dn/dM„ on the reduced mass of the system 
with the Martin potential has been found by the Bohr-Sommcrfeld quantization, so that at the step from 
(bh) to (bq), the density changes less than about 15%. 
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Making the second Borel transform Lj^j2{t) on eq. ( 1U7| ) with account of eq.( p.lU| ), one 



finds the expression for the leptonic constants of the pseudoscalar (6c) states, so that 

^ Ml dF |-^"^np(MJ + Cg(MJ (— G )| , (111) 

where we have used the following property of the Borel operator 

Lr{x) x" exp(-6x) ^ 5^+\t - h) . (112) 



Explicit form for the spectral density and Wilson coefficients can be found in ref. [|76 . 

Expression ( p.ll| ) is in the agreement with the above performed derivation of scaling 
relation (p!7|). 

The numerical effect from the mentioned corrections is considered to be not large (the 
power corrections are of the order of 10%), and the uncertainty, connected to the choice 



of the quark masses, dominates in the error of the fs^ value determination (see eq . (11031) ) 



Thus, we have shown that, in the framework of the QCD sum rules, the most reliable 



estimate of the fs^ value ( |103| ) is coming from the use of the scaling relation ([T7|) for 
the leptonic decay constants of the quarkonia, and this relation agrees very well with the 
results of the potential models. 

3. Decays of Be mesons 

3.1. Life time of Be mesons 

The processes of the Bf. meson decay can be subdivided into three classes a) the h- 
quark decay with the spectator c-quark, b) the c-quark decay with the spectator 6-quark, 
c) the annihilation channel 5+ — >■ l~^i'i{cs,us), I = e, fi, t. 

The total width is summed from three partial widths 

T{B, ^X)= T{b ^X) + r(c ^ X) + r(ann) . (113) 

The simplest estimates with no account for the quark binding inside the B^ meson and 
in the framework of the spectator mechanism of the decay for the first and second cases, 
lead to the expressions 



So, that mfe and rric are chosen for to represent correctly the spectator parts of the total 
widths for the B and D mesons. 

The width of the annihilation channel equals 

r(ann)=5:-^|Hcr/|M,,m^l-^ -C,, (115) 
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The diagrams of the Be meson decays: (a) the c-spectator decay; (b) the 6-spectator 
decay; (c) the annihilation 

where Cj = 1 for the rz/,- channel and Cj = 3|V^sp for the cs channel, and rrii is the mass 
of the most heavy fermion (r or c). 

Note that in the case of the nonleptonic decays, an account of the strong interaction 
results in the multiplicative factor of enhancement to formulae ( p.l4| )- ([TT5|) (see section 
3.2). 

The mentioned widths, calculated with the use of the know values of parameters 
^g, I He I = 0.046, |Ks| = 0.96 etc., are presented in Table [l^. 

Thus, the rough estimate of the life time leads to tb^ ~ (2 — 5) ■ 10~^^ s. So, the fraction 
of the c-quark decay is approximately 50 %, the 6-quark one is 45 %, and the annihilation 



Table 15. 

The widths (10 eV) of the inclusive decays of b- and c-quarks in free and bound 
states (in the Be meson) and the branching ratios {BR in %) of inclusive Be decays 



Decay 


free 


Bt 


BR 


Decay 


free 


Bt 


BR 


mode 


quarks 






mode 


quarks 






6 — >• c + e'^h'e 


62 


62 


4.7 


c — s + e"*" + I^e 


124 


74 


5.6 


b + fi+u^ 


62 


62 


4.7 


c S + /i+ + z/^ 


124 


74 


5.6 


b ^ C + T^Uj. 


14 


14 


1.0 


c ^ s + u + d 


675 


405 


30.5 


b ^ c + d + u 


248 


248 


18.7 


c s + u + s 


33 


20 


1.5 


b ^ c + s + u 


13 


13 


1.0 


c d -\- e^v 


7 


4 


0.3 


6 — > c + s + c 


87 


87 


6.5 


c — i> (i + yU+ + z/^ 


7 


4 


0.3 


b ^ c + d + c 


5 


5 


0.4 


c ^ d + u + d 


39 


23 


1.7 






63 


4.7 


B+ ^c + s 




162 


12.2 


B+ ^c + d 




8 


0.6 


5+ ^ all 




1328 


100 



30 



channel is 5 %. However, these estimates do not take into account a quite strong binding 
of the quarks inside the meson: corresponding corrections to the estimates can reach 
about 40 %. 

Let us consider this effect in the semileptonic modes of decays with the spectator b- 
quark. The final state of such decays generally contains the B^*^ mesons with the more 
small phase space of the lepton pair. 

The effect of the phase space decrease is shown on Figure where the kinematical 
borders of Dalitz plot for the B^ BgC^v decay are compared with the borders for 
the c-quark and calculated at different values of the c-quark mass. As one can see from 
Figure 0, the end-point of the leptonic spectrum is approximately one and the same in 
the different decays 

Ml - Ml 

2Mb, ^ ' 

However, the maximum values of the leptonic pair masses q^ax different. 

One can easily find that the spectator model better describes the semileptonic decay 
D ^ K. In the case of the Be meson decay, the admissible kinematical region is strongly 
reduced. With the account for the phase space reduction in the spectator model, one can 
get H: 

r(5+ ^ XbC+v) ^ 0.71 r(D+ ^ X,e+z/) . (117) 

The effect of the phase space reduction does not sizably appear in the case of decays 
with the spectator c-quark. For such decays, as one can see from Figure ^, the spectator 
model well describes the B meson decays as well as the Be meson decays, and one can 
believe that 

r(5+ ^ X,e+z/) ^ r(fi+ ^ X^e+z/) . (118) 



Another possible way of the estimate is related with the summation of the exclusive 
decays into the channels BgC^u and B*e~^i'. In agreement with the same kinematical 
arguments, their sum is the main fraction of the semileptonic decays |^2[. If one neglects 



the decaying quark momentum inside the Be meson, the admissible region of masses in 
the inclusive semileptonic decay Q Q'eu is varied within the limits 

2 

(mg. + msp)2 < (M^) < mj, + + m,p^ . (119) 

TTlq 



From approximate formula ( |119| ) with the use of the constituent quark masses, one can 
see, that the admissible region in the decay Be — Xe is varied in the limit of 200 
MeV and, hence, the final state is saturated by the lowest states. For the considered case 
{nig = nie = 1.7 GeV, m^/ = = 0.55 GeV and mgp = nib = 5.1 GeV), this region has 
the widths, equal to 340 MeV, that is less than the expected difference of masses between 
the basic state and the first orbital excitation of the (bs) system. 
Thus, one can consider that 



r(5+ ^ Xfee+z/) ^ T{Be ^B, + ev) + T{Be B^ + ev) . (120) 
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The Dalitz diagrams for the semileptonic decays: (1) Be — B^lv, (2) Be — Bgli', (3) 
D ^ K*liy, (4) D Klu, (5) c ^ slv (m^ = 1.7 GeV, m, = 0.55 GeV), (6) c ^ s/i^ 
(rric = 1.5 GeV, = 0.15 GeV); E is the lepton energy, is the square of the lepton 
pair mass 

The results of different quark models for the semileptonic B^. decays (see section 3.2) 



The Dahtz diagrams for the semileptonic decays: (1) B^ iplu, (2) B^ ijjv, (3) 
B — Dlu, (4) B — D*lu, 5) 6 — s- clu, E is the lepton energy, is the square of the 
lepton pair mass 



32 



lead to the following sum of the widths of decays into Bs and B* 

T{B, Bs + eu) + T{Bc B^ + eu) ^ (60±7)10"^^ GeV ^ 0.5 T{D+ X,e+z/) . (121) 

Accounting for the current theoretical uncertainties, one can calculate 

r(fie ^ Xi + e+z/) = (0.6 ± 0.2) r(D+ ^ X,e+z/) . (122) 

For the c-spectator decays, the calculations in quark models and QCD sum rules show, 
that the semileptonic decays are saturated by the transitions into the lowest rjc and J/ip 
states, i.e. 

r(fi+ ^ X^e+z/) ^ r(5+ ^ (r/e + J/^)e+z/) ^ T{B+ X,e+z/) . (123) 
With the account for these factors, the probabilities of inclusive decays are presented 



in Table 0. The widths of the hadronic inclusive decays, which are in details discussed 
in Section 3.3, are also shown. 

The compact sizes of Be meson lead to the large value of the weak decay constant 
if Be ~ 500 MeV), that enforces the role of the annihilation channel into the massive 
fermions c, r. The decays of B^ meson into the light fermions are suppressed due to 
the forbidding over the spirality. Although the use of the effective masses for the u- 
and d-quarks instead of the current masses can increase the width of the annihilation 
channel into ud, the latter will yet be much less than the width into the heavy massive 
fermions. In agreement with eq. (|115|) , the conservative estimates of the annihilation 
decay probabilities are presented in Table |1^. 

Thus, the consideration of three types of processes for the Be meson decay leads to 
the life time estimate 

~ 5 ■ 10-1^ s 

with the following approximate sharing of branching fractions: 37 %, 45 % and 18 % cor- 
respond to the c-spectator mechanism, 6-spectator one and the annihilation, respectively. 

The uncertainty in the estimation of the Be meson life time is generally related to the 
choice of quark masses. The mass of 6-quark rrib = 4.9 GeV is chosen so that one can 
describe the B meson life time in the framework of the spectator mechanism. Note, that 
the differences of the life times for the charged and neutral B mesons are unessential and, 
hence, the given choice of the mass is quite certain. For the D mesons, this is not the 
case, since the life times of and mesons differ twice. 

Nevertheless, there is another way, that is more reliable for the extraction of the c- 
quark mass, this is the consideration of the semileptonic decays of D mesons. Indeed, the 
value rrie = 1.5 GeV in the spectator mechanism well describes the decays — > K^e^u 
and K'e^u, whose widths approximately are equal to each other. Note, at any 

other reasonable choice of rric (from the total widths, say), the error in the Be meson life 
time will be not large, since the summed branching ratio of the Be meson decays due to 
the c-quark decays is about 40 %. 
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Fig. 5. The diagram of the semileptonic decay of Be meson 

3.2. Semileptonic decays of Be mesons 
3.2.1. Quark models 

In the framework of quark models, the semileptonic decays of Be mesons are considered 



in refs. ||30|j3^j3^ . The detailed study of the Be meson decays in the quark model of the 
relativistic oscillator WSB |7T| has been first made in ref . and further in ref . |Q , where 
the quark model ISGW |j70[ has also been used. The covariant description approach, early 
offered for the composed quarkonium model, is developed in ref. |^ . 

Consider the amplitude of the B'^ Mxe^Ue transition with the weak decay of the 
quark 1 into the quark 2 (Figure ^) 



A 



72 



(124) 



where Gp is the Fermi constant, Vu is the element of the Kobayashi-Maskawa matrix. 
The lepton current is determined by the expression 

//. = e(gi)7^(l- 75)^^(92) , (125) 

where qi and q2 are the lepton and neutrino momenta, respectively, (q'l + q'2)^ = t. 
The Hfj^ quantity in eq. (|124|) is the matrix element of the hadronic current 

J^ = V^-A^ = gi7^(l - 75)Q2 . (126) 

The matrix element for the Be meson decay into the pseudoscalar state P can be 
written down in the form 



{Beip)\A^\Pik)) = F+(t)(p + k), + F4t)ip - k) 



(127) 



and for the transition into the vector meson V with the mass My and the polarization A, 
one has 



{Be{p)\MV{k,X)) = -(M + Mv)Ai(t)e(-^) 



M + Ml, 



J^^V) (P + k)^ 
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M + Mv 



Relations ( 127 ), ( 128 ) define the form factors of the Mxe^Uf, transitions, so, 

for the massless leptons, F_ and do not give contributions into matrix element ( |124|) . 
In the Covariant model of the quarkonium (see Appendix I), one can easily find 



FAt) = hnii + mJ^—^p{t), (129) 

2 ^ M 7712 

F-it) = -^(mi-m2 + 2m3p)J^— ep(t), (130) 

2 V -M 777-2 

Here mgp is the mass of the spectator quark (see Figure and the function C,{t) has the 
form 



cj2 + cj|; u^ + ul MMx 



where Mx is the recoil meson mass, lox is the wave function parameter (I.6)-(I.8) for the 
recoil meson. 

tmax = (M-Mx)^ (132) 

^max is the maximal square of the lepton pair mass. 

For the vector state one has My = Mx, and we obtain 



V(t) = l(Af + Afv,)^^-l^v((), (133) 

^ 2 M + Mv V M ma ^ ^ ^ ^ 

^ l(M + M,)(l-?!^)/fi-«,M. (135) 

.43(*) ^ -l(M,M.)(l + 5;^)/fi^?v-(*). (136) 

It is interesting to note, that the exponential form of the form factor dependence on t 
( |131| ) can be quite accurately represented, in the admissible kinematical region, by the 
form, corresponding to the model of the meson dominance 

1 

1 -T/ml 

where are presented in Table [1^. One can see from eqs.( |130| )-( |136| ), that the form 
factors (excepting Ai{t)) are also representable in form ( [137| ), and the degeneration takes 
place 

my = m^a = ~ Tn+ , (138) 



35 



Table 16. The mi. parameters (in GeV) for the (^(t) representation in eq.(137) 







5+ r/cC+z/e 


5+ ^ 5,e+z/e 


B+ ^ B:e+u, 


rrik , GeV 


6.3 


6.45 


1.9 


1.95 



if cjp ^ ujv^ Mp ^ My. 

As for the Ai{t) form factor, it can be represented in the form 



1 



ai + 



A[{0) 



1 - t/mX 



where 



A[{0) 
ai 



my , 

1 + - + 2M{m2 - m^p) 

2 M + Mv 

Ai(o)-a;(o). 



My 1 
M 7712 



ey(o) 



(139) 

(140) 

(141) 
(142) 



The values of the transition form factors at zero mass of the lepton pair are shown in 
Table IT^. The numerical calculations in ref . pO| have been performed for the mass values 



rUb 



4.9 GeV , rrir = 1.6 GeV = 0.5 - 0.55 GeV . 



(143) 



The element of the Kobayashi-Maskawa matrix has been taken equal to Vbc = 0.046. 
The constant in ref. has been varied in the limits 

/^^ = 360 - 570 MeV , (144) 

where the upper limit corresponds to the values, obtained in the nonrelativistic potential 



model p4| , p2| , in the parton model [^, and in the QCD sum rules [|36| , |76| . The lower 
limit corresponds to the value, obtained in the Borel sum rules of QCD |^|7B]. Note, 
that for the — > i^e'^u^ decay, the result weakly depends on the Jb^ choice (3 %). 
One has also supposed 

U = U , (145) 



Table 17. The form factors of the semileptonic Be decays 



mode 




^i(O) 


A'M 


^2(0) 


V{0) 


B+ ^ ^e+u. 




0.73 


0.14 


0.67 


1.31 




0.89 










B+ ^ B,e+v, 


0.61 










Bt - Ble+u, 




0.52 




-2.79 


5.03 
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and one has varied the values 



/bs = 100-llOMeV, (146) 
f^. = 160 - 180 MeV , (147) 



that does not contradict the estimates, made in the QCD sum rules [0 . 

Note, that for the semileptonic Be meson decays Mxe^Ue, where Mx is the 

recoil meson, the explicit covariance of the model allows one to take into account correc- 
tions on the velocity of the Mx meson. As for corrections due to the quark motion inside 
the meson, they are taken into account due to the difference between the constituent and 
current masses of the quark. 

In the ISGW model for the meson state vector, one uses the following nonrelativistic 
expression 

fnig \ _fmq 



\q[—Px + P,s q[^p^-p,rs ). (148) 

where XsP^ is the spin wave function of the quark- ant iquark pair in the state with the 
total spin 5* and the spin projection ms, C*mims is the coupling between the orbital 
momentum L and the total spin S of the system with the total momentum Sx] 4>x{p)LmL 
is the corresponding nonrelativistic wave function, px is the meson momentum, p is the 
relative momentum of quarks. In the considered model, the meson mass is equal to the 
sum of quark masses in the approximation of infinitely narrow wave package, only. 

As the probe functions, the nonrelativistic oscillator wave functions have been chosen 

^ = ^-H-V 

= -l^rexp^ 



The j3 parameters have been determined by the variational principle and the Cornell 
potential 0. 

In the WSB model, the mesons are considered as a relativistic bound state of a quark 



qi and an antiquark q2 in the system of infinitely large momentum ||71|| : 



|P,m,j,j,) = V2(27r)3/2^ /'d=^pid=^P2 5'(p-pi-p2) 

Li^^{pu.x,s,,S2)af{pi)hf{p2m, 

where = (Po,0,0,P), and at P — oo, x = pizjp corresponds to the momentum 
fraction, carried out by the nonspectator quark, pn is the transverse momentum. 



37 



For the orbital part of wave function, the solution of the relativistic oscillator is used 
L^(pt,x) = Nm\lx{l-x) 6^p(~^) 
exp 



2cj2 V 2 2m2 



(149) 



In the both models, the calculation of hadronic matrix elements {Bc{p)\J^\X{k)) comes 
to the calculation of the matrix elements of the quark currents between the quark states 
and the overlapping the corresponding wave functions. 

In the potential models, the bound state of two particles is described by the wave 
function with the argument, being the relative momentum of the particle motion in the 
meson system of the mass centre. However, in the case of the decays with large recoil 
momenta, one can not chose the system, where the both mesons (the initial one and the 
decay product) would be at rest, so that one has a kinematical uncertainty in the form 
factor values. 

For instance, in the ISGW model the form factor dependence on the invariant mass 
of lepton pair t is determined by the function ^(t): 

where j3 and j3i are the parameters of wave functions for the initial and final mesons, mgp 
is the spectator quark mass, M and Mi are the model parameters (the masses of the initial 
and final "mock" -mesons) |j70[. The k parameter in eq. (|150| ) is introduced synthetically 
for the correct description of the electromagnetic form factor of vr meson {k = 0.7). So, 
the authors of ref. related this factor with possible relativistic corrections at large 
recoil momenta. 



Recently in ref. [Q, the model for the description of the heavy quarkonium decays 
has been offered. In this model, the required behaviour of form factors (at k = 0.7) 
is automatical with no introduction of additional parameters. In contrast to the above 
mentioned approaches (the Covariant quark model and ISGW model), the nonrelativistic 
approximation is performed for the hadronic matrix element as a whole, but it is not 
performed separately for the wave functions of the initial and final states, only. At small 
recoil momenta, this formalism practically repeats the ISGW model, but at large mo- 
menta, there are some differences in the structure of the spin part of the wave function 
and the argument of the wave function of the final meson. So, the latter change is the 



most important and it leads to the difference in the form factor dependence on t 



The transition form factors in the ISGW model depend on ^ind Pbs- its 
values, = 0-82 and = O-^l are obtained from the variational principle. Since 
the considered model is the nonrelativistic approximation, the form factors are the most 
accurately predicted at = q'^^^ = {Mb^ — Mx)^ (at the maximal value of the lepton 
pair invariant mass). 

One can calculate the form factors in the region of low values in two different ways: 
by the use of the exponential dependence on as in ISGW or in the pole model of meson 
dominance. The results for the decay widths, calculated in these ways, are presented in 
Table |21|. The additional parameter in the ISGW model is k = 1 ( |150|) . 
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The results, obtained in ref . , is also presented in the same Table. In the constituent 
quark model, the exponential dependence of the form factors can be represented in the 
pole form. As one can see from Table in the ISGW model for the decays, where the 
c-quark is the spectator, the exponential dependence and the pole model give the different 
results. 

In the WSB model, the form factor values at = are predicted dependently of 
the uj parameter (see eq. (|149|) ), that corresponds to the average transverse momentum 
of quarks inside the meson. In ref. |^ the uo values were equal to the average pi values. 



estimated in the ISGW model {ux ~ (3x)- Note, that the uo parameter is external for the 
WSB model. 



The results of the mentioned approaches are presented in Table 
Note, that the relative yield of the pseudoscalar states in respect to the vector states 
is much greater in ref. where r*/r 2 in comparison with r*/r ^ 3-4 in the ISGW 
model. This leads to that, for example, the exclusive decay modes i?^ il){ric)e^Ve 
practically saturate the h — > cev transition. This feature is analogous to the consideration 
of the B D^*^eu decay, that also saturates free fe-quark decay. The decays into the 
excited states and manyparticle modes are suppressed. 

As one can see, these three models for the decays with the spectator &-quark, give the 
close values 

T{Bc Bs + e + u) + T{Bc ^ Bl + e + u) = (60 ± 7) ■ 10~^ eV. 

Note also, that in the case of the heavy quarkonium Be, the application of the non- 
relativistic wave function instead of the wave function of the relativistic oscillator in the 
meson of the WSB model, seems to be more acceptable. This circumstance and uncer- 
tainty in u, maybe, explain, why the WSB model gives the underestimated value for the 
width of the B^ — >• J/ip + e + u decay. 

3.2.2. B+ J/'^{T]c)e+u decay in QCD sum rules 

The most suitable for the registration modes of the Be decays are the semileptonic or 
hadronic transitions with the J/ip particle in the final state. But in the QCD sum rules 
(SR) [^,^,^ and in the quark models, one found different results for both the widths of 
the corresponding decays and the form factors of the transitions, although in the frame- 
work of the separate approach, the calculations, performed in different ways, coincided 



with each other. Recently in ref. [^, one has shown that the existing discrepancy can 
be cancelled by means of account for the higher QCD corrections in SR. 

The widths of the semileptonic B^ decays are defined, in general, by the form factors 
F+, V, Ai and A2 (see eqs.( p.27| ), (|128|) ). Following the notation of ref. [^, redefine the 
form factors ( |127D , ( |128D in the way 

U = F+, F^ = {Mb^ + Mv)A,, 

M V 
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For the calculation of these form factors in the QCD SR, let us consider the three-point 
functions 

^^J,{Pl,P2,q^) = i"^ J dx dy exp{i{p2X - piy)} 

{0\T{c{xh,c{x),V,{0),b{yHc{y)}\0) , (151) 
^lfiPi,P2,q'^) = i"^ J dx dy exp{i{p2X - piy)} 

{0\T{cixMx),J^'\0),biyHciy)}\0) , (152) 

Introduce the Lorentz structures in the correlators 

= n+(pi+p2)M + n-g^, (153) 

^lu = i^V^i.uapP2Pl^ (154) 

Yi% = intg,, + Yi^p^2P'i + ^2P'iPi+'^ip'2P2+^ip'iP2- (155) 

The form factors /+, Fy, and are determined by the amplitudes 11+ , Ily, Hq 
and 11:^ = ^(Hi + 112), respectively. For the amplitudes, one can write down the double 
dispersion relation 

where Q"^ = —q^ > 0. 

The integration region in eq. (|156|) is determined by the condition 

_ . 2gig2 + {si + S2- q^){rnl -ml- Si) 

\y\s,,S2,q^)\^/\mls,,ml) ' ^ ^ 

where A(a;i, 0:2, x^) = {xi + X2 — x^Y — 4x1X2- 

In accordance with the general ideology of the QCD sum rules |TT[, the right hand, 



theoretical, side of eq. ( |156|) can be calculated at large euclidian pi and p^ values by the use 
of the operator product expansion (OPE). The perturbative parts of the corresponding 
spectral densities (the unit operator in OPE) to the one loop approximation are presented 
in Appendix II. Since we consider the systems, composed of the heavy quarks, one can 



neglect the power corrections 36 



Consider the physical part of SR. As has been already mentioned in the consideration 
of the axial constant of B^, there are two approaches. In the first one, one assumes that 
the physical part includes the contribution of the lowest mesons and the continuum, that 
is approximated by the perturbative part of the spectral function from some threshold 
values Sg and Sq [35,^|. The contribution of the higher excitations and the continuum 



is suppressed due to the Borel transformations over two variables {—pf) and (— p^)- The 
numerical results, obtained in such way of ref . |]35| , |36| , are presented below. 

In the second way, one saturates the spectral density by infinite number of narrow 
resonances 13011, so that 
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6{s, - Mhy{s2 - Ml') 

6{s^ - MhM^2 - ) 

ij^i mb + rric 
5{s, - M^yis^ - M;') (158) 

Substituting the expressions for the spectral densities ( |158| )-( p!5^ ) in the dispersion 
relations for correlators (|156|) to the one hand, and their perturbative values to the other 
hand, one gets the corresponding sum rules. 

Applying the procedure, described in Appendix III, for the both sums over the reso- 
nances, one obtains for the form factors under consideration 

2^ _ 8m,(m, + m,) dM^^ dM^ 1 2,2 , 

MMJUZ^^^' ' ' ' Q ) ' (159) 



FoAQ) - -l^mj^^^J^2poAMB.,M^,Q). (160) 

Choosing the k and / values, one can extract the transitions between the given res- 
onances. At k = I = 1, one gets the required form factors for the J/tp{ric)e~^i' 
decays. 

Thus, we use the phenomenological parameters dM^/dk instead of the additional 
parameters such as the continuum thresholds. As has been mentioned, the former is, in 
a sense, the density of the quarkonium states with the given quantum numbers. One can 
quite accurately calculate these factors. The masses of the radial excitations of ip are 
known experimentally [0, and for the and rjc systems, composed of the heavy quarks, 
one can use the predictions of the potential models P- pi]|J5^JF7| -|5B[ . 

The dMk/dk values at A; = 1 for the systems under the consideration, are presented 
in Table 



Let us chose the following values of parameters fs^ = 360 MeV, fr,^ = 330 MeV, pqJTql , 
rub = 4.6 ±0.1 GeV, = 1.4 ±0.05 GeV, gj/^ = 8.1 (from the data on T{J/^ e+e")). 
For the axial constant, we chose 360 MeV instead of 460, to compare the form 
factor values with the results of ref. 0]. The B^ meson mass will be varied from 6.245 
to 6.284 GeV (the data of the different potential models). Note, that at such choice of 
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Table 18. The derivatives dAI^/dk (in GeV) for the lowest states 



Quarkonium 


Be 




Vc 


dMk/dk 


k=l 


0.75 


0.75 


0.76 



the parameters, we do not go out from the integration region ( |157D . In ref. [^| one used 
Mb^ = 6.35 GeV. The form factors values, obtained in ref . p0tp6|,p5| at = 0, are show 



m). As in 



in Table |19|. The deviation from the central values in Table |T9| corresponds to the variation 
of the quark and Be meson masses within the limits, mentioned above (for 
the case of the potential models, the SR predictions agree with each other. 
In ref.pOf, the form factors have the following pole behaviour 



(161) 



pole 



and 



where mpoie = 6.3-6.4 GeV, and 4>i{Q'^) = l + atQ"^. The representation of /+, Fy 
F^ by the form ( |159| )-( P!^nD gives the following Oj values, which are quite low and equal 

considered above, does 
0. The results for 



to -0.025, -0.007, -0.012, -0.02 respectively. The behaviour, 
not practically differ from the ordinary pole behaviour [3D|, where Oj 
the transition widths are presented in Table 

As one can see from Tables ^ and |2^, the results of the Borel SR are, in general, in 
a good agreement with the results of the considered approach, within the model errors. 

the q^- 



The widths, obtained in ref. 35 



are greater, than in ref.[^,0, since in ref. P5| 
dependence of the transition form factors strongly differs from the behaviour, expected 
in the meson dominance model. 

The deviation from the quark models is related, from our point of view, with that in 
the calculations of the transition form factors in the QCD sum rules, one has to account 




Fig. 6. The Coulomb corrections in the semileptonic meson decay 
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Table 19. The form factors of the — > J/tp{r]c)eu transitions at = 



/+(0) 


Fv{0) 
GeV"^ 


Ff(0) 


GeV 


Reference 


0.23 ±0.01 


0.035 ±0.03 


-0.024 ±0.002 


2 ±0.2 


[30| 


0.2 ±0.02 


0.04 ±0.01 


-0.03 ±0.01 


2.5 ±0.3 
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0.55 ±0.1 


0.048 ±0.007 


-0.030 ±0.003 


3.0 ±0.5 
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for as/f -corrections, where v is the relative velocity of the quarks inside the meson. For 
the heavy quarkonia, where the velocity of the quark motion is small, such corrections, 
corresponding to the Coulomb-like interactions (Figure pf), can play an essential role |31 



Indeed, the spectral densities S2, Q^), determining the Be decay form factors, are 
calculated near the threshold si = M^^, S2 = M^^ ^. When the recoil meson momentum 
is small, the calculation of the ladder diagrams in the formalism of the nonrelativistic 

[|17|, Figure leads to the finite renormalization of p, so that 



quantum mechanics (see 



Pi{Si,S2,Q 



max;' 



C Pi{Si, S2, Qmax) 



where the factor C has the form 



C 



B. 



(162) 



(163) 



and •''^'^^(0) are the Coulomb and free wave functions of quarks, so that 



^^(0) 



[1 — exp( 



Ana 



3v 



(164) 



For the two-point quark correlators, determining the decay constants / of the heavy 
quarkonia ip, T, Be, the account for factor (|163|) led to the essential enhancement of /, 
so that one observed the agreement with the experimental data on and fr- Note, that 
the expansion in (|164|) over as/v —>■ exactly leads to the dominant term, appearing in 
account for the one-loop a^-corrections to the two-point correlator of currents. Moreover, 
these corrections have been taken into account in the evaluation of / for the three-point 



Table 20. 

The widths (in units 10 eV) of the semileptonic Be decays in the QCD sum rules 
with no account of a5/i;-corrections 



mode 


m 


Ff 








4.6 


7 


10.5 


Bt 




1.4 


1 


9. 
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correlators, but one did not take into account the loop corrections in the determination 
of the three-point spectral densities. 

One would, in a logics, not to account for the a^-corrections in the / evaluation as 
well as in the p determination, or one would take into account these corrections in both 



cases. As one can see, for example, from ( p.59|) , one can write down 



where the f^^^ and p^^^ values are calculated with no account for the as-corrections, and 
the factors C appear due to the Coulomb-like corrections and defined in eqs. (|163|) , (|164|) . 
It is evident, that 

C 

^1/2^1/2 = 1 • (166) 
Thus, in the determination of the transition form factors, we can use the "bare" / and p 



quantities, calculated in zero approximation over 0:5 |8^, instead of that, say, was done in 
ref. where p^^^ was used without the C factor and the / constants with the account 



for the ag-corrections, i.e. with the factors Cb^ and C^^^^, were instantaneously used. 
As the result, one gets the following values for the form factors /+ and F^^ ||31[ 

/+(0) = 0.85 ±0.15 , F(f = 6.5 ± 1 GeV. 

For the corresponding widths, one has found |^ 

r(5+ ^e+z/) ^ 44 ■ 10"^ eV , T{B+ ri^e^v) ^ 15 ■ 10"^ eV. 

Note, that we have neglected the contributions of the form factors Fy and F^ in the decay 
Be J/ipev. This can result in the overestimation of the widths values up to 10-20 %. 
One can make the agreement between the obtained values of the widths and the results 
of the quark models (see Table ^) within the limits of the theoretical uncertainties of the 
used methods. 

Comparing the results of the QCD SR and the quark models, one can accept as a 
central value of the Be J/ipeu decay width (with the accuracy about 40 %) 

T{Be J/tpeu) ^ 40 ■ lO-'' eV, 

that corresponds to the branching fraction, equal to 3 %. Then the relative probability 
of the three lepton yield in the B^ decays, when two of them reconstruct J/ip, equals 

i?r(i?+^(/+r),//+z/)^8■10"^ 

where /, /' denotes e or p. 

3.2.3. Approximate spin symmetry 

In the bound state, the heavy quark virtualities are much less than their masses, i.e. 
the following kinematical expansion for the quark momentum pq is accessible 

p'^ = mQ-v^ + k'' , (167) 
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Table 21. 

The partial widths (in 10 eV) of the semileptonic Be decays (ISGWl and ISGW2 
are the results of the ISGW model with the exponential dependence of the form 
factors and the pole model, respectively) 



mode 


ISGWl 






ISGW2 g 




WSB @ 












38.5 


53.1 


21.8 


37.3 


34.4 


Bt^ 




10.6 


16.1 


16.5 


20.4 


14.2 


Bt^ 




0.033 


0.12 


0.002 




0.094 


Bt^ 




0.13 


0.32 


0.011 




0.268 


Bt-^ 


^/'(25)e+i/e 










1.45 


Bt^ 












0.727 


Bt^ 




16.4 


17.9 


11.1 


16 ±4 


26.6 


Bt^ 


Ble+Ue 


40.9 


46.3 


43.7 


41 ±6 


44.0 


Bt^ 


Bde^Ue 


1.0 


1.1 


0.5 




2.30 


Bt^ 


B*e+u, 


2.5 


3.0 


2.9 




3.32 



so that 

vk^O, |fc2|<m^. (168) 

Then in the system, where v = (1, 0), the heavy quark hamiltonian in a gluon field of an 
external source has the form 

so that in the limit Aqcd ^ f^q, the spin- flavour symmetry EHQT takes place for 
hadrons with a single heavy quark. 

For the heavy quarkonium, one has purely phenomenologically, that the kinetic energy 
does practically not depend on their flavours, however, the value of the potential energy 
term ^(r) is determined by the average distance between the heavy quarks. This distance 
depends on the quark masses, i.e. the flavours. Therefore, there is no flavour symmetry 
of the wave functions in the heavy quarkonium. However, the magnetic field of the heavy 
quark is determined by its motion velocity (as well as magnetic moment). The quark 
motion is nonrelativistic in the heavy quarkonium, so that 



B ~ Ofv) ~ O 



(—) ■ (170) 



From eqs.( |169| ), (|170D it follows, that the spin-dependent potential in the heavy quarko- 



nium appears in the second order over the inverse heavy quark masses (see Section 2) 



VsD ~ O 



(4) ■ ^'''^ 



Thus, in the leading approximation for the heavy quarkonium, one can neglect the spin- 
dependent forces in comparison with the kinetic energy and the nonrelativistic potential. 
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This means, that in this approximation, the quark spin is decoupled from the interaction 
with the gluons of low virtualities, therefore the masses of the nLj quarkonium states are 
degenerated over J, and these states have the identical wave functions. 

Thus, there is the approximate spin symmetry for the heavy quarks in the heavy 
quarkonium. 

Further, let us consider the matrix element 

M={n''Lj{QQ')\T\h) , (172) 

where F is the operator of quark currents, /i is a state. Then the spin symmetry means 
that the action of spin operators of the heavy quark is factorized, and the matrix element 
M, obtained by the action of the quark Q spin (or by the antiquark Q' spin) 

= \ v^Qa-^ , a-^='- [7- , (173) 

is related with the matrix element M by the equation 

M = {n'Lj{QQ')\S^T\h) = E^ii: {n'' Lj,{QQ')\T\h) , (174) 

where M is the sum of the matrix elements with J', and Cgg, are defined by the rules for 
the spin operator action. 

For the semileptonic ric{ip)l~^i' decays, the spin symmetry is valid in the point of 

zero recoil of rjclip). Indeed, in this case the spectator c-quark and the c-quark, produced 
in the weak decay of the 6-quark, are practically at rest in respect to each other, so binding 
into the state, they interact with low virtualities, characteristic for the heavy quarkonium. 
At a nonzero velocity of the c-quark, it must exchange with the c-quark by a momentum, 
comparable with its mass, for to make the bound state, where their velocities are close. 
Thus, at the nonzero meson recoil, the gluons with high virtualities can shift the heavy 
quark spin, and the spin symmetry does already not take place. 

At the zero recoil of the charmonium vb^ = the covariant amplitude of the 

weak current, the nonzero contributions are given by Ai(t), F±(t) at t = tmax, and the 
heavy quark spin symmetry means, that 

(Mb, + M^J F+ + {Mb^ - M^J F_ = (Mb, + M^) Ai , t = Wx , M,, = M^ . (175) 

Thus, in the approximation of zero spin-dependent splitting of the heavy quarkonium, 
one derives the specific relation for the form factors of the semileptonic exclusive decays 
into the charmonium. 

Note now, that the covariant model, considered above, gives the semileptonic form 
factor values for the Be decay into the charmonium, so that these quantities satisfy the 
symmetry relation ( |1 75| ) . In contrast to the decays of the heavy hadrons with a single 
heavy quark, where the form factor normalization at zero recoil is fixed due to the flavour 
symmetry, the normalization of form factors for the weak semileptonic transitions between 
the heavy quarkonia is determined by the overlapping of their wave functions, which 
depend on the quarkonium model. For the oscillator wave functions in the considered 
potential model, we get 

(Mb, + M^) Ai(t^ax) = PMb2M^ ^(W) , (176) 



46 



where 

In ref.|^^ the factor .^(tmax) was determined in the quarkonium model with the Coulomb 
potential, that is quite rough approximation. 

Note further, that in the semileptonic B^. decay, the lepton pair kinematically has, in 
average, large invariant masses m{l'^v) k, 1.9 GeV, so that the Ai form factor contribution 
dominates, and in accordance with the meson dominance of the t-dependence of the form 
factors, relation ( |175|) , giving Ai{t^.^^^ determines, in a sense, the matrix element of 
the semileptonic ip{'r]c)l^i' decay. This feature can be used for the determination 

of the Be meson mass from the mass spectrum as well as the element \Vbc\ of the 
Kobayashi-Maskawa matrix. 

3.3. Hadronic decays of Be mesons 

Although the semileptonic B^ — > J / tp [i^ {e'^)i' ^{ye) decays can serve as a good trigger 
for the Be registration, the complete Be reconstruction needs a large statistics because of 
the neutrino presence in the decay products. The direct measurement of the Be meson 
mass is possible only in the hadronic exclusive decays. The preliminary estimates of 
some nonleptonic decay widths with the J / ip particle in the final state were made in refs. 
2^,33,^ in the framework of the potential models. The hadronic decays were in details 



considered in refs. [p2| , |34| , p3| . In ref. for the calculations of the transition form factors, 
one used the WSB and ISGW models, mentioned above. In the calculations of the decay 
widths, the phase space reduction was accounted for the c-spectator decays (see Section 



3.1), in contrast to some other calculations [^J81|. In the forthcoming consideration of 



the hadronic decays of the Be meson, we will follow the results of the latter paper. 

The effective four-fermion hamiltonian for the nonleptonic decays of the c- and b- 



quarks has the form [^] 

= ^Vu,,V:jCl{f,)Ol + CUf,)0'L] + h.c. , (178) 
^efr = ^V,,,V*^JC'4^^)0', + Ct{fi)Ot] + h.c. , (179) 

where 

0± = (qialui^ - 75)^/3) (g377''(l - l5)q25){Sa(35jS ± SaS^jfi). 

The factors C±'{fj,) account for the strong corrections to the corresponding four- 
fermion operators because of hard gluons |]3^ , |86| . 

The transition amplitudes must not depend on the subtraction point fi, if one consis- 
tently calculates them in the perturbation theory, i.e. one constructs the corresponding 
functions of the initial and final hadronic states in the perturbation theory, in accordance 
to the operators. 

The problem is complicated, when one deals with the factorization approximation, 
used for the calculation of the matrix elements. In this approximation, one assumes 
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that the current is proportional to a single stable or quasistable hadronic field, and one 
calculates its matrix element between the vacuum and the corresponding asymptotic 
hadronic state, so this procedure gives a value, proportional to the decay constant of 
hadron. After that, the amplitude of the weak decay is factorized and it is completely 
determined by the hadronic matrix element of an other current, that can be calculated 
by use of a model, as in the case of the semileptonic decays. In this approximation, the 
interaction in the final state is neglected. 

Note, that the exact factorization takes place in the leading order of the 1/Nc expansion 
Wl\ . In this approximation, one has to be careful in the choice of the subtraction point. 



since the matrix elements depend on /i. (The dependence of coefficients for the four- 
fermion operators of effective hamiltonian on the subtraction point does not compensated 
by the functions of the initial and final states.) The most suitable choice is /i ~ mc, 
since the radius of the Be meson is determined by the c-quark mass, and the transferred 



momenta in the decays are about nic [B4 



The anomalous dimensions of the 0\ and 01 operators at = rric have the form 
In the leading logarithm approximation at /i > mc, one has [BS] 



as{ml) J \ as{iJt?) 

Clif^) = IWr'- (181) 

at as{ml) = 0.27, as{ml) = 0.19, q;s(M^) = 0.11, one has the values C^(mc) = 0.80 
and (mr) = 1.57. 



At /i > mfe, the anomalous dimensions of the operators are determined by eq. (|180|) , 

but at rric < fi < rUb one finds 



2\ X -3/25 

(183) 



/ ;\^2 ^^ -12/23 ^ . 2^^ -12/25 

Ctip) = ^ (^) . (184) 

The numerical values are C^(mc) = 0.90 and C^lrric) = 1.57. 

For the nonleptonic inclusive spectator decays of Be meson, the enhancement factor 
due to the "dressing" of the four-fermion operators by hard gluons, is equal to 

3lCi^ + C=^l, (185) 

where 3 is the colour factor. For the annihilation decays, the corresponding factor equals 
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The widths of the annihilation and inclusive spectator decays are presented in Table 
r5|. As has been mentioned, the quark masses have the following values rric = 1.5 GeV, 



rrib = 4.9 GeV and = 0.15 GeV, i.e. one takes the choice for the good description 
of both the semileptonic decays of B and D mesons and the total B meson width. The 
enhancement factors ( |185|) , ( |186| ) are calculated in the large Nc limit (this approximation 
gives a good description of B and D meson decays For the fe-spectator decays, 

one accounts for the phase space reduction, in contrast to calculations in ref . . 

The results of calculations for the widths of exclusive decays (here we consider the two- 
particle states) were performed in the models of WSB, ISGW and |]S3| and are presented 



in Tables ^ and The Cabibbo nonsuppressed widths of 6-spectator decays are shown 
in Table g2[ 

The ai and 02 coefficients, accounting for the renormalization of the four-fermion 
operators, are defined in the following way 

iV, + 1 N,-l 
_ iVe + 1 iV,-l 

- ^^^N--^-^N-- ^^^^^ 



In the limit Nr 00 one has 



ai ^ 0.5-(C+ + C_), 

a2 ^ 0.5-{C+-C^). (189) 



The ai and 02 values, used in refs. p^|83| , differ from each other because of a different 



choices for the quark masses and the Aqcd parameter. 

Note, that for the decays with the B mesons in the final state, the contribution of the 
annihilation and "penguin" diagrams is suppressed as 0(sin^°6'c). As one can see from 
Table |2^, the WSB results agree with the ISGW model, and the sum of the two-particle 
decay widths is equal to the total inclusive width of the 6-spectator decay (see Table pIS]) . 



The widths in the model of ref. |8^ are slightly greater. The reason for the deviation 
from the results of two other models can be the fact, that for the 6-spectator decays, 
there are B and Bg mesons in the final state, so that these mesons are relativistic systems 
because of the light quark presence and, hence, the nonrelativistic approximation can 
poorly work. 



Among the c-spectator decays, the widths, presented in Table 23, are given for the 



decays, where the WSB and ISGW models result in close values and one can neglect 



contributions of the annihilation and "penguin" diagrams. As one can see from Table 
the data of ISGW model well agree with the results of model [^ . 

The total inclusive nonleptonic width of the Be meson decay with the J/ ip particle in 
the final state can be obtained from the corresponding width of the semileptonic decay 

r(S, -> J/^X^di,)) = 3 ■ alr{B, ^ J/^ez/)|Kd(s)|' • (190) 

In the large Nc limit, one has 3 ■ af = 4.6 (ai = 1.18) and 

r{Bc ^ J/^X,d(,)) ~ 190 ■ 10-6 eV. 
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Table 22. 

The widths (in units 10 eV) of two-particle hadronic 6-spectator decays {Mb^ = 
6.27 GeV, Mb, = 5.39 GeV, Mb* = 5.45 GeV) 



Decay 


WSB 


ai = 1.23 


ISGW 


ai = 1.23 


m 


ai = 1.12 


mode 




a2 = 0.33 




aa = 0.33 




as = -0.26 




al 31.1 


47.8 


al 44.0 


67.7 


al 58.4 


73.3 


Bt^B, + p^ 


al 12.5 


19.2 


al 20.2 


3.1 


al 44.8 


56.1 


Bt - Bl + 7r+ 


al 25.6 


39.4 


al 34.7 


53.4 


al 51.6 


64.7 


Bt - B: + p+ 


al 115.6 


177.8 


al 152.1 


234 


al 150 


188 


B+ B+ + K'^ 


al 28.2 


3.1 


al 61.4 


6.7 


al 96.5 


4.25 


B+ ^ B+ + 


10.0 


1.1 


al 24.1 


2.6 


al 68.2 


3.01 


B+ B*+ + K° 


al 31.0 


3.4 


al 28.3 


3.1 


al 73.3 


3.23 


B+ B*+ + 


147.1 


16 


163.8 


18 


al 141 


6.23 


B+ ^ B° + TT+ 


0.97 


1.49 


al 1.89 


2.9 


a^ 3.30 


4.14 


B+ -^B^ + p+ 


al 0.94 


1.45 


a2 2.14 


3.3 


al 5.97 


7.48 


B+ 5*0 ^ ^+ 


1.58 


2.42 


al 1.28 


2.0 


a2 2.90 


3.64 


5+ ^ + p+ 


al 8.82 


13.6 


a? 8.86 


12. 


a^ 11.9 


15.0 


5+ ^ 5+ + ttO 


a| 0.48 


0.05 


a| 0.95 


0.1 


al 1.65 


0.074 


5+ ^ 5+ + pO 


a| 0.47 


0.05 


al 1.07 


0.12 


al 2.98 


0.132 


5+ ^ 5+ + u; 


al 0.38 


0.04 


a| 0.87 


0.09 






5+ ^ 5*+ + 71° 


0.79 


0.09 


al 0.64 


0.07 


al 1.45 


0.064 


5+ ^ 5*+ + pO 


al 4.41 


0.48 


4.43 


0.48 


a| 5.96 


0.263 


5+ ^ B*+ + UJ 


a2 3.60 


0.39 


a| 3.53 


0.38 






B+ ^ Bs + K+ 


a? 2.18 


3.35 


al 3.28 


5. 


a2 4.2 


5.27 


B+ B* + K+ 


al 1.71 


2.6 


a? 2.52 


3.9 


al 2.96 


3.72 


B^^B^ + K+ 










a? 0.255 


0.32 


B^^ B^ + K*+ 










al 0.180 


0.226 


Be 5*° + K+ 










a? 0.195 


0.244 


Be 5*0 + K*^ 










a^ 0.374 


0.47 



The branching ratio of the Be decay with the J/ip particle in the final state equals 

BR{Be J/ip + X)^ 0.2 . 



The WSB and ISGW models give close results for the two-particle Be meson decays 
with the B mesons in the final state. Unfortunately, it is complex to detect the Be mesons 
in such decay modes, since one has to reconstruct the B mesons from the products of their 
weak decays. For the Be meson observation and its mass determination, the Be — > J/ ipir'^ 
decay is more preferable |^^. Its branching ratio is equal to 



BR{B+ J/iIjtt+) ^2-10" 
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Table 23. The widths (in units 10 ^ eV) of two-particle hadronic c-spectator decays 



Decay mode 


ISGW 


ai = 1.18 


o 


ai = 1.26 




aj 1.71 


2.63 


aj 2.07 


3.29 




aj 4.04 


6.2 


al 5.48 


8.70 


B+ J/ij + 7r+ 


aj 1.79 


2.75 


al 1.97 


3.14 


B+ J/ij + p+ 


aj 5.07 


7.8 


a2 5.95 


9.45 


B+ -^r], + K+ 


aj 0.127 


0.195 


al 0.161 


0.256 


B+ ^ r/, + K*+ 


a2 0.203 


0.31 


a? 0.286 


0.453 


B+ J/^ + K+ 


aj 0.130 


0.2 


a? 0.152 


0.242 


B+ ^ J/tP + K*+ 


a2 0.263 


0.4 


al 0.324 


0.514 


B, tlj{2S) + 7r+ 






a2 0.251 


0.398 


B,^ij{2S)+p+ 






al 0.71 


1.13 


B^ ij{2S) + K+ 






a? 0.018 


0.029 


Be V^(25) + K*+ 






a? 0.038 


0.060 



The Be meson decays, where the CP violation can be observed, Bf {cc)D^, Be 
Dp{7i) and Be D^Dg are of a special interest. 

The approximate estimates for the decay branching fractions and the asymmetry pa- 
rameter oi CP violation, are obtained in ref. The corresponding results are presented 



in Table ^ The asymmetry parameter A is defined in the following way 



A = — ^-^ ^ — 191) 

A large value of the asymmetry is expected in the Be D^D^ decays with the 
meson, decaying into the eigen CP invariance state. However, the branching ratio of such 
event is too low 

BR{B+ Dl+D°) ^ 10-\ 

The D*~^ meson identification is also complicated. As one can see from Table ^ the 
best mode for the CP violation observation in the Be meson decay can be B^ {cc)D^. 
However, even at the expected statistics of the Be meson yield at the future colliders 
(about 10^-10^^ events), it is difficult to observe such events, if one takes into account the 
branching fractions of the (cc) states and D meson decays. 

It is difficult to estimate the decay widths, but it is worth to mention about the Be 
meson decay modes such as Be 3DX or Be Dgip and Be DK. The Be ip{3S)D 
decay can be of a great interest, when ip{3S) has the decay into the D meson pair. 
However, it is probable, that this decay width, as the width of the decay into three D 
mesons, is small because of a smallness of the phase space, yet. The Be Dgcf) decay 
width can be roughly estimated of the order of 2 % |5^, but it will be very difficult 



to observe the Be meson in such mode because of the complex reconstruction of the Dg 
meson. 



51 



4. Production of Br mesons 



The electromagnetic and hadronic production of Be as the particle with the mixed 
flavour supposes the joint production of the heavy quarks b and c. This explains the low 
value of the Be production cross section in comparison with the production cross sections 
of the particles from the ip and T families. On the other hand, the absence of the B^ 
decays channels into light hadrons due to the strong interactions leads to that all bound 
(be) states basically transformed into the lowest state with the probability close to unit, 
due to the radiative transitions (see Section 2). 

From the theoretical point of view, the production of Be meson, having the small sizes, 
takes place with the virtualities of the order of the heavy quark mass sum. This fact makes 
sure the perturbation theory to be applicable to the processes of the Be production. The 
nonperturbative part, related with the account for the Be wave function, is quite reliably 
calculated in this case. 



4.1. Production of Be mesons in e~^e -annihilation 

The simplest example of the Be production in e"'"e~-annihilation (in the region of Z 
peak) is described by the diagrams on Figure |^. 

The matrix element of the (be) quarkonium production is transformed from the cor- 
responding matrix element T{pb,Pe) for four heavy quark production by the integration 
over the relative momentum of the b- and c-quarks with the weight of the quarkonium 
wave function 

- I ||*«l)nP.P.):^(-« + rn.r-i-P. + ,n.rrf'^^f^^ . (192) 
where M is the meson mass and 

^"^ = i=75"^ rf = ^rAx (193) 



Table 24. The branching ratios {BR) and asymmetries (^4) for the CP violating Be decays 



X 



BR{Bt - X) 



A 



1.5-10-2 
-0.3 ■ 10-2 
0.6-10-2 

10-3 
10-3 
10-3 
10-3 



7]eD*+ 



1.0 
1.2 
0.5 
2.8 
1.6 
3.3 
1.8 
1.6 
0.4 



10-^ 
10-^ 
10-^ 
10-5 
10-5 
10-5 
10-5 
10-6 
10-6 



1.9 
3.0 
1.3 
2.0 



-8.9 ■ 10-3 
-13.8 ■ 10-3 
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Fig. 7. The diagrams of the single Be meson production in e+e -annihilation 



for the pseudoscalar state and the vector one {B^, B*), respectively. The quark momenta 
are determined by the relations 

Ph = jjP + 9, , Pc= jjP - q , (194) 

p-q = 0. (195) 

For the heavy quarkonium one has |q| -C ruf,, rric and eq. (|192|) can be simplified by the 
substitution for T^^{pi,pc) by its value at q = 0. Then 



In papers of refs. plH4^ the total cross sections of the Be and B* mesons and its distri- 
butions over the variable z = 2Eb^I \fs have been obtained. The result of the precise 
numerical calculations in the technique of spiral amplitudes with the Monte Carlo inte- 
gration over the phase space is presented on Figure 8. One can see that this distribution 
is rather hard with the maximum at z^^-^ = M/{M + m^) ~ 0.8. At this value of z^g^^ 
the Be meson and c-quark have zero relative velocity. If one remembers, in the considered 
approximation the relative motion of the c- and 6-quarks inside the B^ meson is absent, 
then one clearly finds, that the maximum in the distribution corresponds to the configu- 
ration, when all quarks move as a whole with one and the same velocity. In this case the 
minimal virtualities of the initial 6-quark p"^ = [nib + 2mc)^ and the gluon k"^ ^ 4m^ are 
realized. At any other z values, these virtualities increase. Note, these speculations are 
correct only for the last two diagrams on Figure ^ when one can neglect the contribution 
of the first and second diagrams, suppressed up to two orders of magnitude in respect to 
the former. In the asymptotic limit s —>■ oo, when one can neglect terms of the order of 
M'^/s and the higher powers of this ratio, choosing the special gauge condition (the axial 
gauge with the four-vector n = (1,0,0, —1) along the direction of the 6-quark motion), 
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one can show that the contribution of the last diagram on Figure ^ survives only. In this 
case the expression for a'^da/dz acquires the sense of the fragmentation function of the 
6-quark into the Be meson, if one chooses the 6-quark production cross section a as the 
normalization factor at the same energy. 

The function of the h ^ Be fragmentation, where Be is the pseudoscalar state, has the 
following form 

8a||^(0)P rzil-zY , 

D{z)i^B. = ^ V n n v^6-181-2r;^ + 

^Iml (1 — (1 — r)z)° 

(21 - 74r + 68r2)z^ - 2(1 - r)(6 - 19r + ISr^);^^ 

+3(l-r)2(l-2r + 2r2)2^) , (197) 

and for the fragmentation into the vector state one has 

Diz)- - 8«.^I^(0)P rz{l-zf 

- 21ml {l-{l-r)zf^^ ^'^'^ 

3(3 - 2r + Ar'^)z^ - 2(1 - r)(4 - r + 2r'^)z^ + 

(l-r)2(3-2r + 2r2)/) , (198) 

where r = mc/{mh + m^). As one can see from Figure H, Di^^^{z) and Di^q*{z) are 
in a good agreement with the results of precise calculations. The b ^ B* process has 
slightly more hard distribution in comparison with the b ^ B^ one. At as = 0.22, 
1^(0) P = fsc = 560 MeV and rric = 1.5 GeV, the corresponding integral 

probabilities are equal to 3.8 ■ 10~^ for b ^ B^ and 5.4 ■ 10~^ for b ^ B*. The probabilities 
of the c-quark fragmentation into B^ are suppressed by two orders of magnitude in respect 
to the values, given above. As the fraction of the bb production, the total number of 
the produced Bc{Bc) mesons with the account of the B*{B*) states and the first radial 
excitations is equal to (at as = 0.22) 

cr(e+e bb) 



Due to the quark-hadron duality, there is the independent way of estimation for this 
ratio. To reach this goal, one has to compare the obtained cross section of the bound be 
state production with the cross section for the production of the colour-singlet (be) pair 
in the process e~^e^ bbec with the low values of invariant mass M^^ 



Ml 



daje+e ^ te) 5e-singl ^ 
dM? 



dMl , (200) 



'be 



where mo = m^ + m^ < Mg^ < Mb + M£, + AM = Mth and AM ^ 0.5-4-1 GeV. Supposing 
mo = 6.1 GeV and Mth=8 GeV as the threshold value, one gets the be system production 
cross section of the order of 7 pb. On the other hand, the sum of the cross sections for 



the production of B^. and its first excitations equals 9.3 pb, as is seen from Table 25 
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Table 25. 



The cross sections (in pb) for the production of the S'-wave states of Be mesons in 
the Z boson peak 



State 






2^50 


21^1 


a{as = 0.22) 


3.14 


4.37 


0.805 


1.078 



The comparison of these two independent estimates states, on the one hand, about 
a good agreement. On the other hand, it means that the contribution of the higher 
excitations is not large, and the total cross section is saturated by the S'-wave levels. 

Recent direct calculations of the cross sections for the P level production [^] confirm 
this conclusion. According to the estimates of this paper, the sum over the cross sections 
for the P-wave levels production is less than 10 % from the sum of the S'-wave level 
contributions. 

In ref. the functions of the heavy quark fragmentation into the heavy polarized 
vector quarkonium have been studied, so, for the longitudinally polarized quarkonium, 
one has found the expression 



8a2|^ 



rz(l 



z 



Bi 



81ml (l-(l-r)^)6 

„2n 2 



;2-2(3-2r)^ + 



(9 - lOr + 16r^)z^ - 2(1 - r)(4 - 5r + Qr^)z^ + 
(1 -r)2(3-6r + 6r2)24) , 



(201) 



that does not depend on the polarization of the fragmentating quark. At r = 1/2, 
expression ( P01| ) coincides with the result, obtained for the heavy quarkonium with the 
hidden flavour (T, 



Fig. 8. The functions of the 6-quark fragmentation into the Be and B* mesons 
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Fragmentation function ( |201| ) agrees with the consideration of the heavy quark frag- 
mentation into the heavy meson (Qg), where in the hmit of infinitely heavy quark, EHQT 
leads to the equal probability production of the vector quarkonium with the arbitrary 
orientation of its spin, i.e. to the absence of the spin alinement and to the ratio of the 
vector and pseudoscalar state yields, equal to \^/P = 3 |l94 |. 



For the heavy quarkonium, the relative yield of the vector and pseudoscalar mesons 
is close to unit, and the spin alinement of the vector state has a notable value. For the 
B* meson, this can be observed in the angle distribution of the B* B^pi decay, which 
composes the total B* width. This distribution has the form 

,_(^l^^^,^e, (202) 



dcos^ V2-^ 

where 9 is the angle between the photon and the B* polarization axis in the system of 
B* rest, and the asymmetry parameter ^ determines the relative yield of the transversally 
polarized 5* state 

i = ■ (203) 

For the integral asymmetry at the small mass of the generated quark, entering the meson, 
r <^ 1, one has 

e = ^ + -|r + 0(0. (204) 
3 lo 

The anisotropy in the B* Bcj decay is numerically equal to 6 %. 

In ref . |^8[ the vector quarkonium spin alinement has been studied versus the transverse 



momentum in respect to the fragmentation axis. One has derived quite bulky analytical 
expressions for the fragmentation functions D^^^* (pt), that hnearly tend to zero at pt — ^ 
and decrease as at pt 00. It is interesting, that the average transverse momentum 
at the fragmentation into the longitudinally polarized vector Be quarkonium is twice 
greater than the average transverse momentum at the fragmentation into the transversally 
polarized B* meson, (pt) ~ 7 GeV. 

The event with the B^ meson has the characteristic signature. The hadron jet from the 
6-quark must be produced in the direction, opposite to the B^ motion. The B^ meson must 
be accompanied by the D meson with the average ratio of the momenta (z^) / {zbJ ~ 0.3 



and the average angle between the momenta about 20° [^] . 

The single Be meson production in e"^e^-annihilation has been also considered in 



refs.ppl 



In ref. [p9] the exclusive production of the B^*^^B^*^ pairs in e^e -annihilation has 



been calculated at low energies, where one can neglect the Z boson contribution. The 
total cross sections of the vector and pseudoscalar states have the form 

3e,f?!!!2 - 1 + A - 3ej2- (1 - "fr '^i> )' , (205) 

V mi / V mi/ mias[4:m2) / 



37 2m6 



a(e+e- ^ iQiQ2)p{QiQ2)v) = " ^'^^r'^^^ /^/^ (l - v'fv' x 
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' m'ias{4:m2 



3^i-3e, ^J^-;^J^ j , (206) 

37 2m§ 



3e, - 3e.^4^^)^[3(l - v^) + (1 + - a)^ + 



'mfa5(4m2) 



a 



2 



^{l-v'){l~3v')], (207) 



where u = ^J^^^^AtvPjl , M = rrii + 



M V eimf as(4mi)// V eimfa5(4m|)/ ' ^ ' 

The relative yield of the B^. meson pairs R = a{B^B^)/a{bb) reaches its maximum at the 
energy ^/s = 14 GeV, where it is equal to i? ~ 10^^. This ratio rapidly decreases with 
the energy growth, where the single Be production becomes to dominate. 

As one can see, the study of the Be meson production in e"'"e~ -annihilation allows one 
to make the analytical researches of the heavy quark dynamics. 

Thus, in the Z° boson pole, where the 6-quark production cross section is large, one 
has to expect of the order of 2 events with the Be production per each thousand bb pairs. 
As it is planned, in the experiments at the LEP accelerator, the number about 2 ■ 10^ Z° 
bosons will be detected. This means, that the total number of Bc{Bc) events has to be 
of the order of 10^. Certainly, the real number of reconstructed events will be less, if one 
takes into account the particular modes of the decay. 

4.2. Hadronic production of B^ mesons 

As has been mentioned, the process of the Be meson production in e"'"e~-annihilation 
at large energies can be reformulated as the process of the b BeiB*) fragmentation, 
appearing with the probability about 10"'^. 

The hadronic Be production turns out to be more complex. First, at hadronic produc- 
tion the region of low partonic energies dominates, so that the asymptotic regime with 
the cross section factorization 

^ ~ D-,_bSz) (209) 

is not yet realized. Second, in the hadron interactions a new type of diagrams appears, so 
we will further label the latter as the recombinational diagrams, for which the factorization 
does not take place. 

The contribution of such diagrams, dominating at low masses of the B'^bc system, 
decreases with the growth of this mass, however, it remains essential even at large masses 
and large transverse momenta. The contribution of such type diagrams into the B^*'^ 
production has been first calculated for the exclusive i?^*-* pair production in the quark- 



antiquark annihilation at low energies |38 



The typical set of the QCD diagrams in the fourth order over is shown on Figure |^. 
Here, as in the case of the Be production in e"'"e~-annihilation, the matrix element of the 
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Fig. 9. 

The diagrams of the single Be meson production in the gluon and quark subprocesses 



(be) quarkonium production is obtained from the corresponding matrix element of four 
quarks production by the integration over the relative momentum of the c- and 6-quarks 
with the weight, determined by the quarkonium wave function. 

At high energies, where the Be production cross sections are accessible for the meson 
observation, the gluon-gluon contribution into the production dominates. 

The energetic spectra of the and B* mesons in the system of mass centre for two 
colliding gluons are shown on Figure |10| at the different values of the total energy ^/s =20, 
40 and 100 GeV. 



The cr{gg Bc{B*)cb) values are presented on Figure |ri| at the several energies of the 
interacting gluons for mi) = 5.1 GeV, rric = 1.5 GeV and as = 0.2. The ratio of the cross 
sections crB*/(JBc is about 3 at the energies 20, 40 and 100 GeV and it is about 2 at the 
energy 1 TeV, when in e"'"e~-annihilation, where the b B^ fragmentation dominates, 
this ratio is ctb'/o'b^ ~ 1-3. 

The variation of the o'B*/crB^ ratio is the consequence of the change in the production 
mechanism. The fragmentatonal component gives the low contribution in comparison 
with the contribution of the recombination diagrams. This can be noted from Figure 
p!0| , where the differential cross sections for the Be and B* meson production, calculated 
by the Monte Carlo integration of the exact expression for the matrix element squared, 
are presented in comparison with the cross section, calculated under the fragmentation 
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Fig. 10. , 

The differential da/dz cross sections for the single production of the Be mesons (a) 

and B* mesons (b) in the gluon annihilation at different values of the total energy 



formulae (p^), (p8|). 



The total cross section of the Bc{B*) mesons is obtained from the partonic one crij{s 



The total cross sections for the single production of Be mesons (empty triangles) 
and B* mesons (solid triangles) in the gluon annihilation in comparison with the 
production cross section (multiplied by the factor 2 • 10~^) of the bb quark pairs (solid 
line) 



59 



by the convolution with the functions of the parton distributions in the initial hadrons 

s l—s/s 

ds J dx ^ , / , , 4s\V2 

4(r?i6+mc)2 -1+s/s 



^tot{s)= J ^ J J E/:(^i)/^(^2)a.,(5), x* = (x^ + ^) . (210) 



The cross sections, calculated with the account of the known parameterizations for fal{x) 
|9^ , are presented in Table p6 . 



The energy 40 GeV is close to the c.m.s. energy for the carrying out the fixed target 
experiments at the HERA accelerator. At y/s = 1.8 TeV we present the cross section 
of the Be production in pp-coUisions at Tevatron, and, finally, the energy y/s = 16 TeV 
corresponds to the conditions of the pp-experiment at the future LHC collider. The 
energetic dependence of the cross section, summed over the particle and antiparticle Be 



production, is shown on Figure 12 



From the values, presented in Table it follows that at y/s = 40 GeV, the summed 
cross section o"sum for the meson production is about 10~^ of the total cross section of 
the bb production, so this makes the Be study practically to be impossible in this experi- 
ment. One has to note, that in this case we can not restrict ourselves by the gg —>■ B^cb 
contribution and we have taken into account the contribution of the qq B^cb process. 

The experiments at Tevatron and LHC, where cTsum/cbb is about 10~^, will give the 
real possibility for the observation of hadronic B^ production. Therefore, at the energies 
of these two facilities, we present the most interesting distributions of the cross sections 
for the l^So and 1^5*1 states production (note, that as our calculations show, the cross 
section at the energies under consideration is completely determined by the gluon-gluon 
interaction, since the quark-quark contribution is suppressed by two orders of magnitude, 

10-2). 

The distributions for the 1^5*0 pseudoscalar and l^Si vector mesons are shown on 



Figures ^ and |I4| at the energy of the interacting hadrons 1.8 TeV. 

The distributions da/dx (see Figure [T^^) show, that we deal with the central B^ 
production, where the complete cross section is collected in the interval from —0.3 to 0.3. 
The average transverse momentum of Be is about 6 GeV, and from the distribution over 
the angle between the directions of the Be and c-quark motions, one can conclude that 



c-quark generally moves in the direction close to the Be one . 

One has to note, that the considered diagrams of the QCD perturbation theory are 
the diagrams of the fourth order over as- This results in the strong dependence of the 



The cross sections (in nb) of hadronic production of the BeiB*) mesons (the standard 
deviation in the last digit is shown in the brackets) 









2'So 


2^S, 


atoti^O GeV) ■ 10^ 


1.63(2) 


9.5(2) 


0.13(1) 


0.75(2) 


atotilOO GeV) ■ 10^ 


7.8(2) 


36(1) 


1.1(2) 


5.2(2) 


atot{l.8 TeV) 


13.3(8) 


53(3) 


2.7(2) 


10.4(5) 


atotilQ TeV) ■ 10-2 


1.96(8) 


7.6(2) 


0.43(2) 


1.66(8) 



60 



The total cross sections (in nb) for the single production of Be mesons (empty circles) 
in pp-interactions at different energies and the cross sections (in mkb) for the beauty 
particle production (solid circles) 



cross section on the particular as choice. The latter must be determined by the typical 
virtuality in the production process. The analysis shows, that this virtuality is large 
in the contributions, decreasing faster that 1/s, only. In the remaining contributions, 
including the fragmentational one, it is not large and about Arricmi,. Under this reason 
the as = 0.2 value, chosen as the strong coupling constant, is the most reasonable at this 
scale. The use of the running coupling constant as{s), for example, leads to the decrease 
of the Bc{B*) production cross section by about 7 times. The pessimistic estimates of the 
Bc{B*) production are presented in ref. [Q, where, as one can see, the as{s) value has 



been used. 

At low energies of hadron collisions, the quark-antiquark annihilation with the Be pro- 
duction dominates in respect to the gluonic one, since, in this case, the latter has a much 
lower luminosity, that decreases also with the growth of the total energy of the partonic 
subprocess. At the low energies of the quark-antiquark annihilation, the exclusive B^B^ 
pair production can be essential. The total cross sections of the vector and pseudoscalar 
Be meson production due to the quark-antiquark annihilation have the form 

a(l-,l-) = ^4^'v/l^(1.3 + 1.4A + 0.3A2), (211) 

<^(0-.0-) = ^ ^ A'VT^ (1 - , (213) 
wherefrom one can see, that the vector state production dominates. In eqs.( pil| )-( |213|) 
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Fig. 13. 



The differential da/dpt cross sections for the single production of Be and B* mesons 
in pp interactions at the energy 1.8 TeV 



we have introduced the notations 



A = 4MVs, 

ITlljUlc 
^ = 

TJib + rric 



The differential cross sections for the single production of Bf, and B* mesons in pp- 
interactions at the energy 1.8 TeV; (a) da/dy, where y is the particle rapidity, (b) 
da/dx, where x = lEjs^l'^ 
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Table 27. 

The total cross sections (in units 10 pb) for the pair production of Be mesons due 
to the quark-antiquark annihilation in pp{pp) interactions at low energies 



GeV 




a(l-,0-) 


a(0-,0-) 


30 


0.9 (0.08) 


0.24 (0.022) 


0.006 (0.0004) 


40 


5.8 (0.94) 


1.6 (0.25) 


0.054 (0.007) 


50 


15.8 (3.5) 


4.3 (0.95) 


0.18 (0.034) 



The numerical estimates of the total cross sections for the B^. production in pp inter- 



actions are presented in Table 27 



Summing up the consideration of the hadronic production, one can draw the conclu- 
sions. 

1. The mechanism of the hadronic Bc{B*) production strongly differs from the pro- 
duction in e"'"e~-annihilation. 

2. The relative fraction of the fragmentation contribution is low even in the region of 
large transverse momenta. 

3. The vector state production is enforced in respect to e^e^^-annihilation. 

Thus, the hadronic Be production requires the analysis for the large number of dia- 
grams and its detailed study opens the possibility for the investigation of effects of the 
heavy quark dynamics in the higher orders of the QCD perturbation theory. As for the 
Be yield at the real physical facilities, it is quite high, but the registration of the Be events 
is essentially determined by the detector acceptance (cuts of the transverse momenta of 
particles, characteristics of the vertex detector an so on). 

4.3. Be meson production in uN-, ep- and 77- collisions 

In the previous sections we have considered the Be meson production in the processes, 
where one has the maximal current statistics for the production of hadrons with heavy 
quarks, i.e. at Fermilab and LEP colliders. In the present section we consider the estimates 
for the Be meson production in the processes of the deep inelastic scattering of neutrino 
and electrons by nucleons and in 77-interactions at future facilities. 

4.3.1. Be production in z/A^-interactions 

The diagrams of the neutrino-production of the Be mesons on quarks and gluons are 



show on Figure 15 



Note, that for the Be production in the neutrino collisions with gluons, the suppression 
of the partonic subprocess cross section by the factor | HcP in comparison with the partonic 
subprocess of the Be neutrino-production on light quarks, is compensated by the more 
high luminosity of the gluonic subprocess in comparison with the quark one, so that 
the both mechanisms of the Be production in the neutrino-nucleon scattering give the 
comparable contributions, and (j{vN —>■ BeX) ^ 10"''^ cm^ at the neutrino energy E^, ^ 
500-1000 GeV in the laboratory system. 
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Fig. 15. 

The diagrams of the Be meson production in the processes of the neutrino scattering 
on gluons (a) and quarks (b) 

After the integration over the valent parton d distribution, the c-quark production in 
the W*~^d — >• c process, suppressed as sin^ 9c, has the value, comparable with the c-quark 
production in the W*~^s — > c process, since the strange quark "sea" is suppressed in 
respect to both the valent quark distribution and the "sea" of the more light d-quark. 

The estimates of the Be meson production cross sections, calculated on the basis of 



the diagrams on Figure ^ agree with the estimates, obtained in the model of the vector 
meson dominance (Figure 0) and in the model of the soft gluon emission of the (be) 
pair, that in the colour-singlet state and with the low invariant mass M{bc) < Mb + M^,, 
transforms, in accordance with the quark-hadron duality, into the [be) bound state, which 
radiatively decays into the basic ISq state in a cascade way with the probability, equal to 
1. 

As a result, one can reliably state that the total cross section for the Be meson pro- 
duction in z/A^-collisions is of the order of 10~^ from the total cross section of the uN- 
scattering, so that, at a characteristic statistics about 10^ events in neutrino experiments, 
one can expect only several events with the B^ meson production. 

4.3.2. Production of B^ mesons in ep-scattering 

In contrast to z/iV-scattering, in ep-collisions in addition to the processes of the weak 
charged current exchange, the main contribution into the B^ meson production will give 
the processes with the virtual 7-quanta exchange (Figure ^7^. 

The exact calculation of the diagrams on Figure ^ is not yet performed at present. 
However, one can think that the estimate, made in the Monte Carlo simulation system for 



hadron production HERWIG [49|, is quite reliable, since the HERWIG parameters have 



Fig. 16. The diagrams of the B* meson production in the model of vector meson dominance 
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been chosen to get correct values for total hadronic cross sections of the charmed and 
beauty particles production, being in an agreement with the experimental values. More- 
over, the HERWIG estimates of the production cross sections in e~^e~ and hadronic 
interactions agree with the values, obtained in the exact calculation of the diagrams in 
the QCD perturbation theory. 

Thus, in accordance with the estimates in the HERWIG system, one can expect about 
10^ events per year with the production at the HERA facilities. This Be yield is 
comparably close to that of at LEP. However, the extraction of Be events at HERA is 
complicated by the presence of a hadronic background, that is essentially lower at LEP. 



4.3.3. Photonic production of Be mesons 

At present, the future 77-colliders with the high luminosity (~ 10^^ cm~^s~^) are 
intensively discussed. In this section we calculate the cross section of single Be production 
at energies y/s about 30 GeV in accordance with the diagrams, shown on Figure |T8|. The 
calculation technique coincides with that of described in the section on the hadronic 
production of Be. 

The total cross sections of the Be and B* production are presented in Table where 
one takes as ~ 0.2. One can see that near threshold the pseudoscalar state production is 
suppressed in comparison with the production of the vector one, so at ^/s = 15 GeV one 
has (Jb'/o'Bc ~ 55. Such behavior of the (Jb*/o'Bc ratio has been noted in [6], where the 
strong suppression of the pseudoscalar meson pair production in respect to the vector one 
takes place in the quark-antiquark annihilation. At large energies of the initial photons 
this ratio decreases and becomes aB*/cFBc ~ 4. The inclusive cross sections aB^ and aB* 
have the maximum at ^/s = 20 — 30 GeV and with the s growth they fall like the total 
cross section for the heavy quarks a^^i production. 

The distributions cr~^dcr/d2 over the variable z = 2|p|/y^, with p being the meson 
momentum, are shown on Figures |TU| and ^ for the Be and B* mesons. As follows from 



these figures the scaling in these distributions is broken: with the energy growth the shift 
onto the low z values takes place. Note, the analogous picture has been observed in the 
gluonic production of Be mesons. 

Note, that the detailed consideration shows, that in the matrix element of the 
77 — > bbcc process and hence in the 77 Bebc matrix element one can distinguish 
three groups of contributions, which are separately gauge invariant under both the gluon 
field transformation and the photon one. The first group of contributions is composed 
of the diagrams when the quark production is independent (we will label these diagrams 



Fig. 17. The diagrams of the Be meson production in parton process of 7*(7-scattering 



65 



Fig. 18. The types of diagrams in the photonic production of B^. mesons 



as the recombination diagrams), the second group consist of the diagrams, where the 
(cc) pair is produced from the 6-quark hne (we will mark these diagrams as the 6-quark 
fragmentation diagrams, their contribution will be denoted as cr^"^^"'^), the third group 
contains the diagrams with the (bb) pair production from the c-quark line, so that they 
are c-fragmentation diagrams with the corresponding contributions denoted as a'^~f^"-9_ 



In refs. [^Pli p3| the assumption was offered that the 6- fragmentation contribution has 
to dominate at large values of the Be transverse momentum, independently of the type of 
the process. So the approximate equation has to be valid: 

d_^K!Z= j ^(^).^^-^d,, (214) 
dPt J dkt z z 

2Pt/v^ 

where daqq/dkt is the differential cross section for the production of the fragmenting q- 
quark in the Born approximation, kt is its transversal momentum, and Dg^Bci^) is the 
function of the q —>■ Be + X fragmentation. 

Remind that in the e"'"e--annihilation the 6-quark fragmentation dominates and the 
c-quark fragmentation contribution is suppressed by two orders of magnitude. In the 77- 
interactions, the c-quark fragmentation contribution is enlarged due to the quark charge 
ratio {Qc/QbY — 16 and, therefore we can not neglect it (as one does in e"'"e--annihilation). 
Note further, that the c-quark fragmentation contribution and the 6-quark fragmentation 
one are related to each other by the simple permutation of the quark masses and charges 
(mc ^ mfc and Qc ^ Qh) ( |214| ). 

The distributions dcr/dpt, da'^-i^'^^/dpt and / dpt at 100 GeV for the B^ and 

B* meson production are shown on Figures ^ and The distributions, predicted in 
accordance with eq.( pi4| ) for the 6-fragmentation (curve 1) and c-fragmentation (curve 
2) are also shown. One can see that, as in the hadronic production, the contribution of 
the recombinational type diagram is essential at any reasonable values of the transverse 
momentum of Be meson and it can not be neglected, when one calculates the cross sections 



Table 28. The cross section (in pb) of the photonic production of B^. {B*) 



v^, GeV 


15 


20 


40 


100 




5.1 ■ 10-3 


3.8- 10-2 


6.7-10-2 


2.5 - 10-2 




2.8 ■ 10-1 


6.0-10-1 


4.0 - 10-1 


1.1 - 10-1 
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Fig. 19. 



The cross section distributions, normalized to the unit, over z for the Be meson pro- 
duction at different energies 



even at the large transverse momenta. One can see from the figure, that for the h- 
fragmentation contribution over the greater than about 30 GeV, the fragmentational 
mechanism gives correct predictions. Thus, in its maximum at the energy 20-30 GeV 
the total cross section, including the 5* and corresponding antiparticle production, is 
about 1 pb. This corresponds to 10^ Sc, produced at the 77-collider with luminosity of 



Fig. 20. 

The cross section distributions, normalized to the unit, over z for the -B* meson 
production at different energies 
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The da /Apt, Aa"-^'"'^ / dpt and da"-^'''^ / dpt distributions over the transverse momen- 
tum for the invariant contributions into the cross section of the Be meson production 
at 100 GeV. The curvers 1 and 2 correspond to the prediction of the fragmentational 
mechanism (217) for the 6-quark (1) and c-quark (2) 



The da/dpu dcr^'^'^^s/dpt and da^-^""^^ / dpt distributions over the transverse momen- 
tum for the invariant contributions into the cross section of the B* meson production 
at 100 GeV. The curvers 1 and 2 correspond to the prediction of the fragmentational 
mechanism (217) for the 6-quark (1) and c-quark (2) 
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lO'^'Vjm^^s^^. At large energies, the cross section falls like the bb pair production one. 
The Be production mechanism is close to that of in the gluon-gluon interactions, and it 
does not come to the simple 6-quark fragmentation. 

5. Conclusion 

The discovery and study of the family of (be) heavy quarkonium with the open charm 
and beauty will allow one significantly to specify the notion of the dynamics of the heavy 
quark interactions and the parameters of the Standard Model of elementary particles 
(such values as the b- and c-quark masses, the coupling of the b- and c-quarks - jVbd 
etc.). The present review tends to the aim of a creation of a theoretical basis for the 
object-directed experimental search and study of the {be) heavy quarkonium family. 

Summing the considered problems, one can note the following. 

We have shown that below the threshold of the (6c) system decay into the BD meson 
pair, there are 16 narrow states of the Be meson family, whose masses can be reliably 
calculated in the framework of the nonrelativistic potential models of the heavy quarko- 
nia. The flavour independence of the QCD-motivated potentials in the region of average 
distances between the quarks in the (bb), (cc) and (be) systems and their scaling prop- 
erties allow one to find the regularity of the spectra for the levels, nonsplitted by the 
spin-dependent forces: in the leading approximation the state density of the system does 
not depend on the heavy quark flavours, i.e. the distances between the nL-levels of the 
heavy quarkonium do not depend on the heavy quark flavours. 

We have described the spin-dependent splittings of the (be) system levels, i.e. the 
splittings, appearing in the second order over the inverse heavy quark masses, Vsd ^ 
Oll/mbrric), with account of the variation of the effective Coulomb coupling constant of 
the quarks (the interaction is due to relativistic corrections, coming from the one gluon 
exchange) . 

The approaches, developed to describe emission by the heavy quarks, have been applied 
to the description of the radiative transitions in the (6c) family, whose states have no 
electromagnetic or gluonic channels of annihilation. The last fact means that, due to 
the cascade processes with the emission of photons and pion pairs, the higher excitations 
decay into the lightest pseudoscalar Be meson, decaying in the weak way. Therefore, 
the excited states of the (be) system have the widths, essentially less (by two orders of 
magnitude) than those in the charmonium and bottomonium systems. 

As for the value of the leptonic decay constant Jb^ , it can be the most reliably estimated 
from the scaling relation for the leptonic constants of the heavy quarkonia, due to the 
relation, obtained in the framework of the QCD sum rules in the specific scheme. In the 
other schemes of the QCD sum rules, it is necessary to do an interpolation of the scheme 
parameters (the hadronic continuum threshold and the number of the spectral density 
moment or the Borel parameter) into the region of the (be) system, so this procedure 
leads to the essential uncertainties. The /b^ estimate from the scaling relation agrees with 
the results of the potential models, whose accuracy for the leptonic constants is notably 
lower. The value of /^^ essentially determines the decay widths and the production cross 
sections of the Be mesons. 
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The theoretical consideration of semileptonic Be decays shows, that the results of the 
potential quark models agree with the predictions of the QCD sum rules, if one accounts 
for the Coulomb-like a5/f -corrections. In this case, the approximate spin symmetry 
in the sector of heavy quarks allows one to derive the relations for the form factors of 
semileptonic Be decays in the rest point of the recoil meson. 

The Be meson production allows in some cases the description on the level of analytical 
expressions, such as the universal functions of the heavy quark fragmentation into the 
heavy quarkonium. The fragment at ional mechanism dominates in the B^ production in 
the e^e~-annihilation at high energies (in the peak of Z boson) and it can be studied at 
the LEP facilities. 

The hadronic production of Be is basically determined by the processes of the b- and 
c-quark recombination, since the partonic subprocesses have the most large luminosity in 
the region of low invariant masses of the produced system {bbcc). The Be meson yield in 
respect to the production of the beauty hadrons is of the order^ of 10~^. 

Modes of Be —>■ ipX decays with the characteristic signature of the J/ip particle have 
the quite large probability 

BR{B^ ^ 'ipX) ^ 0.2 . 

Therefore, the Be particle search can start from the separation of the events, containing 
the J/ijj particle, whose production vertex is beyond the primary intersection point. The 
selected set of the events will, of course, contain the background from decays of ordinary 
heavy- light B mesons {bu, bd, bs), since the probability of the B J/ipKX decay is 
about 1 %, and the heavy-light B meson yield is three orders of magnitude greater than the 
Be one. The background separation requires the down-cut over the effective mass of the 
J/ijjX system, where X denotes the charged particles, having tracks from the J/ip vertex. 
The most preferable channel for the Be extraction is that of the B'^ — > ipl'^vi decay, since 
Be is the only heavy particle with the three lepton vertex of the decay ipl'^ V^l'^l^ . 
The probability of this channel is equal to 

BR{B+ tpl^ui) ^ 8 % , l = e, fi, T . 

At a quite large statisticsf the events with the decay B^ — > ipl'^i' can give the possibility 
for the determination of the Be mass value under the ip^ mass spectrum or the missed 
transverse momentum of neutrino in respect to the direction of the Be motion (see Figure 



231) . The necessary condition for the such measurement is a quite high separation of 
charged hadrons and leptons. 

The straightforward measurement of the Be mass can be made in the mode of the 
B^ J / Tpn^ decay, having the branching ratio, equal to 



BR{Bt ^7r+) ^ 0.2 



^In the present review we do not consider in details the Be production in the neutrino-nucleon inter- 
actions, where one can expect only several events with the Be production per year, since the coupling 
constant of the h- and c-quarks is low , so that these processes have no practical significance for the 
experimental search of Be- 

•^The CDF facility with the vertex detector at the Tevatron FNAL has, in this sense, a preferable 
position. 
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The distribution over the invariant masses of the ipl (a) and tplumis (b) systems in 
the — > decay, where Umis is the neutrino with the momentum, equal to the 

missed transverse momentum in respect to the direction of the Be meson motion 

The mode of the —>■ J /ipTi^Ti^Ti^ decay, where three vr mesons can compose the ai 
meson, also is of the interest. This mode must have a significantly greater probability, 
than the Be J/ip'u decay. 

Since the B^. production at the colliding e^e~ beams takes, as mentioned, the frag- 
mentational character, in general (see Figure it must be accompanied by the D meson 
presence in the same jet, where the Be candidate is being observed. Such signature of the 
event would turn out to give a large advantage for the B^ meson search at e^e~ colliders 
in respect to the search at hadron colliders, where the recombinational mechanism dom- 
inates in the B^ meson production at the accessible energies in the nearest future (see 
Figure H). However, one must take into account the possibility of that the probability of 
the 6-quark fragmentation production of the free cc-quark pair is one order of magnitude 
greater than the probability of the fragmentation into the B^ meson and the single free 
c-quark. This means, that with account for the branching ratios for the B and B^ decays 
into J/ipX, the events with the B^ decay and the single D meson will appear only two 
times more often, than the decay of the heavy-light B meson into J/ ipX with the instan- 
taneous production of two D mesons in the same jet. It is not clear, whether one can 
quite effectively separate these two processes at the present of the vertex detectors, i.e. 
whether one can not lose the vertex of the second D meson. 

It is evident, that the progress in the experimental study of the Be meson and general 
physics of the heavy quarks will be mainly related with the development of the vertex 
detectors, so that the latter would give the possibility of a reliable observation of several 
heavy quarks instantaneously (to search the cascade decays, for example). However, since 
at the present statistics of LEP and Fermilab, several dozens of the Be meson production 
events must be observed, one can think, that the practical registration of Be will be 
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realized in the nearest future. 
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6. Appendices 

I. Covariant quark model 

Consider the general statements of the covariant description of the composed quarko- 
nium model. 

By definition, the energy fraction, carrying out by the quark i in the {QQ') meson, is 
its constituent mass rrii, so that 

M = m + m', (I.l) 

where M is the meson mass, m and m' are the fixed values. For the four-momenta, one 
has 

where P is the meson momentum, q is the relative momentum of quarks inside the meson. 
For the quark propagator, one has 

S{k) = {k^Y + m) D{k) . (1.3) 

The constituent quark has, in fact, the fixed energy, so that in the D{k) function, only 
the imaginary part gives the contribution. In the meson rest frame, one has 

'^m D{k) = — 6{\ko\ - m) . (1.4) 
m 

Eq. ([T^) with account for eq. ([r2|) can be rewritten in the covariant form 

ttM 

'^mD(k) = 5(Pq) . (1.5) 

m 

The quark-meson vertex can be represented as 

L,gM = m^v'ik') D-\k) D-\k') xiP; q) , (1.6) 

where v and v' are the quark spinors, the D{k) function is defined in eq.(p3D, F is the 
spinor matrix, determining the quantum numbers of meson. 

The nonrelativistic description of the meson means, that the form factor is determined 
by the expression 

x(P;g)=27r5(Pg)0(g2). (1.7) 
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In the following, we suppose 

0(g2)=iV exp(^) . (1.8) 

The choice ( |1.8D reflects the typical form of the S-wave functions of the charmonium 
and bottomonium, and it allows one to perform the analytical calculation of the semilep- 
tonic decay widths for Be meson. 

Let us define the decay constants / for the pseudoscalar and vector mesons 

{0\J,,{x)\P{q)) = ifpq^exp{tqx} , (1.9) 
{0\Ux)\V{q,\)) = z/v My exp{zgx} , (I.IO) 

where A is the vector meson polarization, and the quark currents are 

■Jbi^ix) = Qix)'j5'j^Q'{x) , (1.11) 
J,{x) = Qix)j,Q'ix), (1.12) 

In the nonrelativistic potential model, one has 

fp^fv = f, (1.13) 

so that 

/ = 2^vl>(0), (1.14) 

where ^1/(0) is the quarkonium wave function at origin. The oscillator function, resulting 
in eq. ([TBI) , has the form 

= f^Y" . (1.15) 



,27r/ "^V 4 

Condition ([I.14| ) means that the normalization constant in eq. (p8|) equals 

iV ^ A f . (1.16) 

mm J 

Thus, for the quark-meson form factor, one finds 

x{P;q) = 27r5{Pq) ^ ^ exp ^ , (1.17) 

mm' J \uj^ / 

where u is determined by eqs. (|I.14D and ( [1.15|) , so that the only free parameter of the 
model is the constant /. For the ip particle, can be, for example, related with the 
width of the ip ^ e~^e~ decay 

r(^ - e+e-) = ^ al^ el ^ , (1.18) 

where Cc = 2/3 is the -quark electric charge. From eq.( [1.18| ), the experimental value of 
the leptonic width [|l5l gives 

= 410 ± 15 MeV . (1.19) 
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As for the fs^ and /^^ values, these constants are determined theoretically in the 
framework of the QCD sum rules and in the potential models. 

Note, that the stated model of composed quarkonium gives, for instance, the exact 
formula of the nonrelativistic Ml-transition for the electromagnetic decay of the vector 
state into the pseudoscalar one V P-y 

1 fi 

YiV ^P7) = Y ^? ' (1-20) 
where is the 7-quantum energy, and the magnetic moment /i equals 



1 , / e e' 

where e and e' are the quark electric charges in the units of the electron charge. 



(1.21) 



II. Spectral densities for three-particle functions 

The spectral densities for the three-particle functions are determined in the following 
way p6| 



p+(Si,S2,Q^ 



3 (k,^ 



A, 



A;[m3(m3 — mi) + 



m3(m3 - ma)] - [2(siA2 + S2A1) - u{Ai + A2)] 



u 



irin 



+ mim2 — m2ms — mim^ 



(11.1) 



Pv{si,S2,Q'^ 



{(2S1A2 - MAi)(m3 - m2) + (2S2A1 - UA2) 



^3/2 

(m3 — mi) + m^k 



(11.2) 



Po{si,S2,Q^) 



^1/2 



(mi - m2) 



, 2 Ai 
m2[m^ - — 



ml + t(siA2 + S2AI - nAiAa) 

rh 

A2 



mi I mg 



mg 



mo 



-(Ai + A2 - m) + mim2 



(11.3) 



P+(Sl,S2,Q^ 



^3/2 



mi 



2S2A1 - UA2 + 4A1A2 + 2A^ 



+ 



mimg 



4s2 - 2u 



+ m2 



2S1A2 - mAi 



2(3s2Ai + S1A2) - 



m(3A2 + Ai) + k + 4A2A1 + 2A^ + m^(4s2 - 2m) + 

6 r 

-(mi - m3) 4S1S2A1A2 - m(2s2AiA2 + siA^ + ssA^) + 
2s2{siAl + S2Aj) 



(11.4) 
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where 

k = (si + S2 + Q^)^-4siS2 , 

U = Si + S2 + , 

Ai = si — ml + ml , 
A2 = S2 - ml + ml . 

In the Be r]c{J/ip)eh> decays, one has mi = and m2 = m^ = m^ for the masses. 



III. QCD sum rule scheme for three-point correlators 

Let us consider the sum rules for the form factor 

1 



00 



Be 



'i 2 T\jj 2 

Be fj ^'^^vc fii(rfl\ 



mb + mc 2m, 



dsids2. 



(III.l) 



Applying the Borel operators L^-^ (— Pi) and (— P2)? defined in Section 2, to eq.( pi.l|) , 
one derives the following sum rules 



2{mf, + mc)m, 



J dsids2P+(si, S2, Q^) exp(-siri - S2T2) (111.2) 



Introduce the notation 



= E fiMCniQ') exp(-M;V2) 



(III.3) 



and transform the left hand side of eq.( [111.2| ) with the use of the formula by Euler- 
MacLaurin pO 



00 r An 

i=l , ,1. Be 



i=l 

n=k—l 



(111.4) 



n=0 



Acting by Lr'{{M^J^) to eq. (|111.2|) and accounting for eq. ([111.4| ), one gets 



T.fiMCf^'mexp{-M^\2) = 
_ 2mc{mb + mc) dM^^ 2 



I p{Ml\s2,Q')eM-S2T2) (111.5) 
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Making the analogous procedure for the sum of the rjl resonances, one obtains 

fki((j2. _ 8meK + me) dM|^ dM^ 1 2 , 2 , 

^+ ^ - M^M^/^j; dfc dZ {2nyP^^^^^ ' ' ^ ^ ^"'-'^ 

Here we have used the property of the Borel operator 

L^(x)(x"exp(-te)) ^ S^+\t - b) . 

It is not complex to generalize this procedure for the remaining form factors. 
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